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Abstract: 

We perform an extensive study of the scale dependence of flavor-singlet contributions to 
the structure function g2{x,Q'^) in polarized deep-inelastic scattering. We find that the 
mixing between quark-antiquark-gluon and three-gluon twist-3 operators only involves 
the three-gluon operator with the lowest anomalous dimension and is weak in other 
cases. This means, effectively, that only those three-gluon operators with the lowest 
anomalous dimension for each moment are important, and allows to formulate a simple 
two-component parton-like description of g2{x,Q'^) in analogy with the conventional de- 
scription of twist-2 parton distributions. The similar simplification was observed earlier 
for the nonsinglet distributions, although the reason is in our case different. 



* Permanent address: Department of Theoretical Physics, Sankt-Petersburg State University, St.- 
Petersburg, Russia 

t Unite Mixte de Recherche du CNRS (UMR 8627) 



Contents 

T Introduction 



2 The Operator Product Expansion 




2.1 Twist-2 




2.2 Twist-3 




2.3 Transverse spin parton densities: why and why not 





3 Hamiltonian approach to the three-particle evolution equations 

3.1 Coefficient functions of local operators 

3.2 Conformal symmetry 

O Helicity basis 

3.4 Conformal basis 

p. 4.1 Ceneral construction 

p. 4. 2 Quark-gluon conformal basis 



3.5 QCD evolution kernels 



p.5.1 


Diagonal evolution kernels 




|3.5.2 


Utt-diagonal evolution kernels 


3.6 Expansion in conformal operators 




|3.6.1 


Quark coefficient function 




t5.6.2 


Cluon coefficient functions . 



4 Results 








4.1 Diagonal evolution kernels 




4.1.1 The highest quark- ant iquark-gluon state 


4.1.2 The lowest quark-antiquark-gluon state 





4.1.3 'i'he lowest three-gluon state 



4.2 'I'he quark-gluon mixing 

4.3 Reduction of the mixing matrix: The two-channel DCLAP equation . . . 



5 Summary and conclusions 



A 


Appendix: 


Calculation of the anomalous dimensions 










B 


Appendix: 


Kacah symbols for the SL[2, R) group 












C 


Appendix: 


Evolution kernels in the conformal basis 



References 



1 Introduction 



Twist-three parton distributions in the nucleon are attracting constant interest as unique 
probes of quark-gluon correlations in hadrons. Quantitative studies of twist-three effects 
are becoming possible with the increasing precision of experimental data at SLAG and 
RHIC, and can constitute an important part of the future spin physics program on high- 
luminosity accelerators like ELFE, eRHIC, etc. 

The structure function g2{x,Q^) of polarized deep inelastic scattering received most 
attention in the past. The experimental studies at SLAG |jl|, Q have confirmed theo- 
retical expectations about the shape of (72(3;, Q^) and provided first evidence on the most 
interesting twist-3 contribution. On the theoretical side, a lot of effort was invested to 
understand the physical interpretation of twist-3 distributions (see e.g. ||, |] for the 
review of various aspects) and their scale dependence []7], |, ^, |TT], |T2|. One- loop 
corrections to the coefficient functions have been calculated []T3|, |T1|, |T^]. 

In spite of the significant progress that has been achieved, understanding of the scale 
dependence of g2{x,Q'^) still poses an outstanding theoretical problem. The difficulty 
is due to the well-known fact |^, ^ that the structure function g2{x,Q'^) presents 
by itself only one special projection of a more general three-particle quark-antiquark- 
gluon correlation function in the nucleon that depends, generally, on two variables - the 
momentum fractions carried by partons. 

In deep inelastic scattering with a transversely polarized target, only this special 
projection can be measured. On the other hand, the scale dependence of the parent 
quark- ant iquark-gluon correlation function involves the "full" function in a notrivial way 

and the knowledge of one particular projection g2{x, Ql) at a given value of Ql does 
not allow to predict g2{x, Q"^) at different momentum transfers: a DGLAP-type evolution 
equation for (72 (x, Q"^) in QGD does not exist or, at least, is not warranted. The reason is 
simply that inclusive measurements do not provide complete information on the relevant 
three-particle parton correlation function. 

From the phenomenological point of view this conclusion is not satisfactory since it 
would mean that one cannot relate results of the measurements of g2{x, Q"^) at different 
values of to one another without model assumptions. The theoretical challenge is, 
therefore, to find out whether the complicated pattern of quark-gluon correlations can 
be reduced to a few effective degrees of freedom. In this case one will be able to find a 
meaningful approximation to the scale dependence that introduces a minimum amount 
of nonperturbative parameters. 



Such an approximation is known [10| for the flavor-nonsiglet (NS) contribution to 
the structure function. To explain this result, it is convenient to use the language of 
the Operator Product Expansion (OPE), see Sect. 2 for more details. The statement of 
the OPE is that odd moments n = 3,5, .. . of the structure function g2{x,Q'^) can be 
expanded in contributions of multiplicatively renormalized local quark-antiquark-gluon 
operators^ 



{{Ot,(M, (1.1) 



We neglect the Wandzura-Wilczek twist-2 contributions throughout this paper. 



where C^_3 are the coefficient functions and {{0^_^{fi))) are reduced matrix elements 
normahzed at the scale /i; 7^.3 are the corresponding anomalous dimensions that we 
assume are ordered with k: 7°_3 < 7^.3 < . . . < 7^13 for each n, and b = ll/3Nc — 2/3nf. 
Note that the number of contributing operators rises linearly with n in the r.h.s. of (|1.1|) . 
This should be compared with the familiar case of leading twist-2 distributions. There, a 
single operator (fiavor-nonsiglet) exists for each moment. A measurement of the moment 
of g2{x, Q"^) cannot separate between contributions of different operators and is, therefore, 
not sufficient to predict the scale dependence. 

The situation simplifies drastically, however, in the large limit. It turns out [|l^ 
that the tree-level coefficient functions, C^_3, of all operators other than the one with 
the lowest anomalous dimension for each n are suppressed by powers of l/N^ so that 
to this accuracy one can approximate the sum in ( |1.1| ) by the ffist term /c = 0. The 
corresponding anomalous dimension 7^^° can be calculated analytically and result can 
be reformulated as a DGLAP-type evolution equation 



4Cp 



+ S{1 



c + ' 



2Cp, (1.2) 



1 - z_ 

where Cp = {N^ — 1) / {2Nc) . Here, we have also included the l/N^ corrections calculated 



m 



12 



The present paper is devoted to the extension of this analysis to the fiavor-singlet 
sector in which case twist-3 composite local three-gluon operators have to be included. 
Coefficient functions of three-gluon operators vanish at tree level, so that gluon contri- 
butions appear entirely through the evolution. The number of independent three-gluon 
operators is, roughly speaking, half of the number of quark-gluon operators and, similarly, 
is rising with the moment n. The subject of this work is to find out whether the whole 
set of gluon operators contributes significantly, or one can reduce the gluon contribution 
to a certain single degree of freedom. We find that such a reduction is indeed possible 
and formulate a two-channel DGLAP-type evolution equation for the structure function 
g2{x,Q'^) that presents our main result. 

In physical terms, the approximation constructed in this paper corresponds to the 
introduction of quark and gluon transverse spin parton distributions which are identified 
with the particular components of quark-antiquark-gluon and three-gluon parton correla- 
tion functions that possess the lowest anomalous dimension. We will find that, ffist, the 
structure function g2{x,Q'^) is dominated by contributions of these two distributions at 
large scales (including the 0{as) gluon contribution calculated in |]15|), and, second, the 
'leakage' of transverse spin to genuine three-particle degrees of freedom at lower scales is 
small due to the specific pattern of the QCD evolution. We would like to emphasize that 
importance of this result is not so much in the possibility to calculate the scale depen- 
dence, but in the identification of important transverse spin degrees of freedom that are 
preserved by QCD interactions. 

The outline of the paper is as follows. Sect. 2 is mainly introductory. We introduce 
necessary notation and present a summary of the existing calculations of the coefficient 
functions in the operator product expansion. Sect. 3 is devoted to the general formal- 



ism of the renormalization of twist-3 gluon operators. We emphasize importance of the 
conformal symmetry and introduce a convenient framework that allows to treat renor- 
malization as a quantum mechanical problem with hermitian Hamiltonian. This section 
is necessarily rather technical and a reader who is only interested in the applications may 
prefer to skip this discussion and go over directly to Sect. 4 where we collect our results. 
The main result of this paper is the generalization of the DGLAP evolution equation 
( p..2| ) to the flavor-singlet channel; it is repeated in Sec. 5 which also contains a summary 
and conclusions. 

In Appendix A we present a detailed calculation of the relevant anomalous dimensions. 
In Appendix B we collect necessary formulae for the Racah 6j-symbols of the SL{2, R) 
group. The conformal basis representation of the QCD evolution kernels is discussed in 
Appendix C. 



2 The Operator Product Expansion 

The hadronic tensor which appears in the description of deep inelastic scattering of 
polarized leptons on polarized nucleons, involves two structure functions^ 

Wl,f = ^e,,^pq^ {s^9i{xB. Q') + ^ - ^/l g2{xB, Q')] (2.1) 

^ P-q I L p-Q J J 

and is related to the antisymmetric part of the Fourier-transform of the T-product of two 
electromagnetic currents: 

W^.^^ = -ImT,!:^) 



iTlt^ = '-j d'xe^'^^p,s\T{j,{x/2)U-x/2) ~ Ux/2)j,i-x/2)}\p,s) . (2.2) 

The light-cone expansion of ( f^.2D at — goes in terms of nonlocal light-cone operators 
of increasing twist, schematically 

^T{jf,{x)ju{-x) -ju{x)jf,{-x)} ^■ = ° 

{ [C * Of^r-' + [C * O^r-'' + (higher twists)} , (2.3) 

where C * Ojs stands for the product (convolution) of the coefficient functions and oper- 
ators of the corresponding twist. This expression is explicitly f/(l)-gauge invariant, i.e. 
d/dx^T{j^{x)j^{-x)} = 0. 

2.1 Twist-2 

For the sake of completeness and in order to facilitate the comparison to twist-3, we 
collect here the relevant portion of the results for the leading- twist. 



^We define the nucleon spin vector as — u(p, s)7^75m(p, s) where u{p, s) is the nucleon spinor 
u{p, s)u{p, s) = 2M, so that = -iM"^. 



Retaining the tree-level quark contribution and the leading gluon correction that 
starts at order 0{as) one obtains ||15 



[C * O^r-' = ^ J2el fdu I \q{ux) ^i,q{-ux) + {x ^ -x)] (2.4) 

q=u,d,s,... ^ "i^i 

_l — ^in(_2;2^^-j _l_ (m Inu + u{l — u))Tt {Gx^{ux)Gx^{—ux)}^ , 

where G^^ = Xr^G^^t"" and the subscript [. . .]^2 indicates the normalization scale of the 
operator. Here and in what follows it is implied that the gauge invariance of nonlocal 
light-cone operators is restored by including the gauge factors Pexp {i J dz^A>'{z). Note 
simplicity of the answer ( |2.4| ): the entire gluon contribution can be eliminated by choosing 
the proper scale of the quark operator = l/(— x^e^'^^). This property is lost in the 
momentum space since after the Fourier transformation contributions of different light- 
cone separations get mixed. 

Going over to the matrix elements, one introduces the quark and gluon helicity dis- 
tributions 

iiim{-x)\p,s) = (sx) j d^e^'^P''Aq{^,fi^), 

{p,s\Tt{g,o,{x)G,^{-x)]\p,s) = '-{sx){px) j\^^^^iAg{i,^?). (2.5) 

In the first case, positive (negative) ^ correspond to the contribution of quarks (anti- 
quarks) Ag(xB) = q^xs) — q^xs), Ag(— x^) = Ag(xs) = q^i^xs) —q^ixs), respectively. 
For gluons, Ag{^) = Ag{-^). 

Ignoring the gluon contribution for a moment, making a Fourier transformation of 
( p.4| ), taking imaginary part and comparing with the definition of structure functions in 
TTD one obtains 

giixB, Q') = IY1 [^^(^^' = + ^9i-XB, f^' = Q')] , (2.6) 



2 



where xb = Q^/{'2pq) is the Bjorken variable, and twist-2 contribution to g2{xB,Q'^) 

[92ixB,Q')r-' = gV'^ixB.Q') = -giixB,Q') + f -giiy,Q') • (2.7) 

Equation (|2.6|) tells that the structure function gi{x, Q"^) is a measure of the quark helicity 
distribution in the nucleon, as well known. Equation (|2.7| ) is the familiar Wandzura- 
Wilczek relation |T^. We would like to stress that despite the fact that the relation in 
( [^.6D is affected by higher order perturbative radiative corrections, the Wandzura-Wilczek 
relation between the structure functions, Eq. (|2.7|), is exact to all orders of perturbation 
theory and is only modified by twist-3 contributions to g2{xB,Q^) that are subject of 
this paper. The reason for this is that the relation ( |2.7| ) follows from particular (and 
unique) form of the Lorentz structure for the antisymmetric part of the T-product in 
( p.3|) which in turn is dictated by the U{1) gauge invariance. Although the coefficient 



function in front of the twist-2 operator in ( p.3| ) has a nontrivial perturbative expansion, 
it affects both structure functions gi{x, Q^) and g2{x, Q"^) simultaneusly and in the same 
way. Hence the Wandzura-Wilczek relation holds true. 

For most of the subsequent discussion it will be convenient to go over to the moments 
space: 



Jo 



for any function /. In particular, restoring the gluon contribution we get |]15| 



[92{n,Q')] 



2.2 Twist-3 



2Mtw-2 



q=u,d,s,... ^ 

a, n — 1 



27r n{n + 1 
n — 1 



-Ag{n, /i^ 



In 



MS 



n 



gi{n,Q' 



(2.9) 



The twist-3 contribution to the T-product in ( |2.3| ) is more complicated. One obtains the 
following expression |]l5|p| 



[C * 0/3]*""^ 



— y 



q=u,d,s, 



1 PU 

du dv ■l{u + v) Sp{u, V, —u) + (m — f )S'^(— -u, f , u) 



+u 



2 



^{u + v)Op{v, u, —u) + 2uOi3{u, V, —u) + (m — f )0/3(m, — m, v) 



(2.10) 



4tt r r ~ ~ 

H — ^ (ln(-x^/i|jg) + 27^; + l) (mu + i-u^ + mIum) vOp{y,u,-u) + uOi3{u,v,-u) 

-vdp{u,-U,v) --^M^(m + 2) 0/3(m, -M,t;) + 0/3(M,t;, -m) + 0/3(t;,M, -m) | , 

where u = 1 — u and we have introduced the C-even nonlocal quark-gluon operator 

Sp{u,v, -u) = S^{u,v, -u) + Sp{-u,v,u) , 

S'J(a,6,c) = ]^q{ax)[igG p^{hx) ± gG i3,^{hx)-i^] ^q{cx) (2.11) 
and the nonlocal three-gluon operator 



Op{u,v,w) = ^r'^Gl^{ux)Gip{vx)GUwx) . 



(2.12) 



^To avoid misunderstanding, note that this result does not include the ©(a^) correction to the 
coefficient function of quark-antiquark-gluon operators. This correction has been calculated recently in 



is], using a different operator basis 



n 



The relatively complicated expression in the last two lines in ( |2.1CI| ) reflects quark-gluon 
mixing and reproduces the corresponding term in the renormalization group equation for 
the twist-3 operator Sf^{u, v, —u) |P, [1^ 

,2 /-n 



[Sf3iu,v,-u)]^2 = [Sp{u,v,-u)]^2 - -^In^ / ds dt{2u 



,-3 



X |[2u(s -t)+ A{u - s){t + s)]Op{s, v,t) - 2\u\{s - t)[0^(s,t, v) - Op{v, s, t)]|.(2.13) 

Similar to the twist-2 contribution in Eq. ( p.4|) this contribution can be eliminated by 
appropriate choice of the normalization scale of the quark operator. Finally, the expres- 
sion in the second hne in (|2.10|) is not affected by the scale choice and defines a 'genuine' 
twist-3 gluon coefficient function. 

Nucleon matrix elements of the nonlocal operators in ( p.l2| ) define parton 

correlation functions 

{p, s\S^{u, V, -u) \p, s) = 2i{px) [s^{px) - p,{sx)]j^ Vi e'v-li^-^i^-i^D^ii,, 6, 6) 

(2.14) 

and 

(p, s\d^{u, V, -u)\p, s) = -2{pxf [s^{px) - p,{sx)]j'v^e'P-^^^'^+^'^-^-^-^Dg{^,, ^3) 

(2.15) 

where the integration measure is given by 



-1 



(2.16) 



-1 



The quark correlation functions D'^l^i) have the following symmetry property: 

^J(6,6,e3) = (Djn-^s, -^2,-^1). (2.17) 
They are in general complex functions, but the imaginary parts do not contribute to 



the structure functions and can be omitted [T8|. In turn, the gluon correlation function 
Dg[^i) is real and antisymmetric to the interchange of the first and the third argument: 

^9(6,6,^3) = -^^,(^3, 6,^1)- (2.18) 

Substituting (|2.1CI| ) into ( |2.3| ) and taking the Fourier transform one obtains the moments 
of the structure function for positive odd n 



[g2{n,Q 



2Mtw-3 



q=u,d,! 



n 



(2.19) 



+ 



a. 



MS' 



An n+1 



fx- 



ip (n) -7ij- 1 



MS 



Here -Dg(^j) is the distribution function corresponding to the C— even combination of 
quark-gluon operators (|2.11|) 



{p, s\S^{u, V, -u) \p, s) = 4z{px) [s,{px) - p^isx)] j e^^'^K^"+«^^-«^"]D,(ei, 6, ^3) 

(2.20) 



so that 

^,(6,6,^3) = ReD+(ei,6,e3) =ReD;(-e3, -6, -eO, (2.21) 

the quark coefficient function is defined as Q 

f) tn-l _ \n-l 

niii.i^) = = , if (2.22) 

and the gluon coefficient functions can be expressed in terms of $^ as 

m^.) = ['^'^1(^1,6) + K-ii-^i - ^3,6)] + (ei -6) , (2.23) 

^n^e.) = (l + ^;^) *'n(e.) + [*Ll(ei, - 6) + K-li-C3, + 6)] . 

Exphcit expressions for $^ and $^ for the lowest moments n = 3, 5, 7 can be found in 



Equivalently, one may choose to expand moments of the structure function in contri- 
butions of local operators, for example 

[G,]% = ^gr'^Glc. (d -xf Gl^ (B -xr-^'' Gl^ . (2.24) 
The reduced matrix elements ((...)) of local operators 

{P,s\[Si]%\p,s) = 2{zpxf^' [s,{px)-p,{sx)] {{[S^U , 

{p,s\[G,]%\p,s) = 2{zpxf^' [s,ipx)-p,isx)] {{[G]%)) (2.25) 
are equal to moments of the three-particle distributions!] 

fvaiit'D^iii) , mN)) = j\e.it'-'mi) ■ (2.26) 



The symmetry relations (CT) and ( CT) imply that (([5^]^))* = (-l)^(([5^]j^~'^)) and 
(([G]^)) = — {{[G]^~^~^f) . Therefore, the number of independent quark and gluon matrix 
elements contributiong to a given n = {N + 3)th moment is equal to £g = + 1 and 
ig = [A^/2], respectively. Finally, we define the reduced matrix elements corresponding 
to the C-even quark-gluon operator ( |2.11| ) as 



%)) = Re {{[S^U = {-If Re {{[S-]^~')) = fva\it'm^) ■ (2. 



27) 



^In difference to Ref. we prefer to define the quark contribution in terms of ReD^(^i) instead of 
RcDq"(^i). Because of tliis, the quark coefficient function in (2.22) differs from the one given in by 
the replacement 1^3 <-)■ —^1. The gluon coefficients are symmetric under this transformation and are not 
affected. 

^The distinction between 'plus' and 'minus' distributions is delicate since it is affected by the con- 
vention used to define the 75 matrix. We use 75 = «7o7i7273 = — i7°7^7^7'^ and e"i23 _ —eQi2z ~ 1 p9[ |. 
A sign change in the definition of the 75 matrix results in the replacement ^ S~ and, for matrix 
elements (([S"]^)) (— l)^(([S']j^~'^)). The coefficient functions in the OPE are changed accordingly, so 
that the results for the physical observables remain intact. 



n 



2.3 Transverse spin part on densities: why and why not 



In order to pursue a parton model-like interpretation, one can introduce transverse spin 
distributions as the specific projections of general quark-antiquark-gluon and three-gluon 
operators. Their definition is suggested by the explicit form of the contribution of these 
operators to the operator product expansion (|2.10|) 



dv{p,s\ [{u + v)S^{u,v,-u) + {u- v)S^{u,v,-u)] \p, s) 



—I 



\sx] 



(px) 



(2.28) 



and 



dv{p,s\ {u + v)0^{v,u, —u) + 2uO^{u,v, —u) + {u — v)0f3{u, —u,v) \p,s) 

= 2 [s,{px)-p,{sx)]j\e'^^^^^ e A^?T(e,/^') 



(2.29) 



so that 



I'd^c-'Aqrio = 4^W*^(e.)^.(e.), 

I'diC-'^griO = 2 I've $^(6) D,{C,) . 



(2.30) 



The functions ^qT^XB.Q"^) and AgT{xB,Q'^) describe the momentum fraction distribu- 
tion of the transverse spin of the proton and have the same support property as the parton 



distribution in ( |2.5|) . Note that Agxi^) = AgT{—^)- However, in contrast with (|2.5| ), 
they do not have any probabilistic interpretation but rather can be expressed through 
the more general three parton correlation functions -Dg(^i, ^2, ^3) and Dg{Ci,^2,^3) inte- 
grating out the dependence on one momentum fraction. The explicit expressions for the 
lowest moments ( p.30| ) look as follows 



4 / VxDg{xi) = 4: 



'oh 



f dx (Aqrix) + Aqri-x)) = 0, 
Jo 

/ dx x'^{AqT{x) + AqT{—x)) 
Jo 

J dx x'^{AqT{x) + Agr(— x)) ~ ^ J (-^i ~ 2x1X3 + 3xl)Dg{xi] 

= 4[{{[S]l))-2{{[S]l)) + 3{{[S]l))] 



(2.31) 
(2.32) 
(2.33) 



for the quark distribution and 



/ dx AgT{x) = / dx x^Agxix) = 0, 
Jo Jo 

1 .1 

dxx^Agrix) = 10 VxxiDg{xi) = -10(([G]^)) 
Jo 



(2.34) 
(2.35) 



for the gluon distribution. 

At tree level, neglecting the (9 (as)— correction to ( |2.19|) , one obtains using (|2.30|) 

™(X5, Q2)]*--3 = \Y.^lf ^[AgHy) + AgT(-y)] (2.36) 

that looks very similar to the leading twist expression ( ^.6p and ( p^.T) ). The two con- 
tributions in the square brackets can be interpreted as the contributions of quarks and 
antiquarks, respectively. Following the analogy with the leading twist, it is convenient to 
introduce the combinations of definite signature 

Aq$iy, Q') = Aqriy, Q') ± Aqri-y, Q') . (2.37) 

The distribution Aq^{y, Q"^) corresponds to the even signature and can be obtained by the 
analytic continuation from even moments N of the OPE. The distribution Aq^{y,Q'^) 
can be obtained by the analytic continuation from odd moments N and defines the 
valence quark contribution. The gluon contribution enters into ( p.l9| ) through 0{as) 
corrections and, according to ( p.30|) , its 'genuine' twist-3 part is parameterized by the 



gluon distribution AqtIx). This suggests that similar to the leading twist expressions 
( ^^ and ( p.7|) , the twist-3 structure function g^~'^{x) can be described in terms of the 
quark and gluon distributions, Aq^{x) and Aqt^x), respectively. 

A deficiency of this interpretation is, however, that it does not go through beyond 
the leading order. This is seen explicitly on the gluon contribution in Eqs. ( p.l9[ ) and 
( p. 23 ): The coefficient function Q'jfl^i) that is responsible for the mixing with quark- 



gluon operators does not coincide with $^(^i) and, therefore, this mixing brings in gluon 
contributions that are not expressed entirely in terms of Ag-ri^) defined in (|2.3CI|) . An- 
other reason is that the scale dependence of the distributions introduced in ( ^.28] ), (|23 



involves the full three-particle functions Dq and Dg in a nontrivial way and, again, brings 
in additional degrees of freedom. 

Aim of this paper is to analyse the effects of QCD evolution in some detail. We 
will find that although the above mentioned difficulties do exist, their numerical impact 
is likely to be minimal. We will then be able to write an approximate effective two- 
channel evolution equation involving the two distributions Ag^(?/, /x^) and Ag^iy, fJ^"^) in 
full analogy with the flavor-singlet DGLAP evolution equations in the leading twist. 

3 Hamiltonian approach to the three-particle evolu- 
tion equations 

Choosing = one can eliminate large logarithmic corrections to the gluon coef- 
flcient function in ( p.l9[ ). To the leading logarithmic accuracy {LL) and retaining the 



flavor-singlet (S) contribution to the structure function g2{x) we write 

IT- J -I 



where 



) = — 5Z , Dl{i,) = Du{ii) + Da{ii) + Ds{^,) + . . . . (3.2) 



Uf 

■' q=u,d. 



The gluon distribution is not present explicitly (to this accuracy) but reappears through 
the evolution of the quark distribution to lower scales. To see how this happens, expand 
in contributions of flavor-sing let local operators defined in ( ^1271) . Using 

and (I^ID one finds 



2 ^ 



g\\n, Q') = {el) - 5^(-l)^-'=(iV - k + 1){{S'AQ'))) ■ (3-3) 



n 

k=0 



The scale dependence of the reduced matrix elements is described by the system of 
coupled evolution equations 

Q"^2((GNiQ'))) = {{V^'U'iiS'kiQ'))) + mXm'{{G^'{Q')))) (3.4) 

with [V^'^] being the known matrices of anomalous dimensions k ,k' = 0, ...,N and 
m , m' = 0, [A^/2] — 1. (Here = n — 3 is the number of derivatives in the quark-gluon 
operator, [. . .] stands for an integer part.) Solving these equations one defines [3A^/2] + 1 
linear combinations of the matrix elements 

{{On,.))= J2 C^kiN) {{[S]%)) + J2 CUN){{m), a = 0,...,[3iV/2] 

0<k<N 0<m<[N/2]-l 

(3.5) 

that are renormalized multiplicatively and obtain the moments of the structure function 
in the standard form ( p..l|) . The corresponding anomalous dimensions 7^ can be found 
by diagonalizing the full matrix of the anomalous dimensions entering 



{CI CD [V^ - 7^1lJ = , = J . (3.6) 

The left eigenstates of the mixing matrix define the vector of the coefficient functions 
entering (13). 

For lowest values of the moments the mixing matrix Vat looks as follows. For = 
the matrix consists of only one element 

while for = 2 it has the following form 



_AA7-l--3-J__|_ln 3Ar_|_23J__|_2 287 Ar 37 1 , 1 
20 ^^c -r 20 Afc ~'~ 5 5 -r 20 TVe ~'~ 5 60 60 iVc ~'~ 5 



120 c 40 c 120 c 3 60 ' 

(3.8) 



-1 -1 



Going through an exphcit calculation of ( |3.5| ) and ( p.6| ) and putting Nc = Uf = 3 one 
finds for the two lowest moments ( p.3|) 



l{el)-'9^\3,Q')=L'r'^/\{S',)), 

lieyg^'-i^, Q') = L^'/' [oA15{{S',)) - 2.558{{Sl)) + 2.776((5°)) + 0.966((G°)) 

+ L^2Aro.34o((^2)) - 0.152((S2')) - 0.261((5°)) - 0.114((G°)) 

+ L-^2/4o.i34((^2)) + 0.m{{S^)) + 0.404((S'°)) - 0.909((G'°)) 

+ L^'/''\o.ni{{Sl)) + 0.2U{{Sl)) + 0.080((S°)) + 0.057((G°)) (3.9) 

where L = as{Q'^)/as{fi'^), all reduced matrix elements on the r.h.s. are normalized at 
the scale /i^ and the flavor-singlet anomalous dimensions are equal to 

7° = 10.5556 , 7° = 10.7393 , -fl = 13.5155 , 7^ = 17.6794 , 7^ = 18.1714 . 

(3.10) 

If /i^ = Q2^ then L = 1 and the coefficients in front of ((^|)), ((^2^)), ((^2°)), ((^2)) coincide 
with their tree-level values 1, -2, 3 and 0, respectively, as expected from ( |3.3| ). Note that 
the largest coefficients occur in the contribution of the operator with the lowest anomalous 
dimension, that is similar to flavor- nonsinglet case 0], and the most important correction 



is apparently associated with the operator with the second-largest anomalous dimension. 
Aim of this work is to explain this structure and understand how it extends for arbitrary 
moment n. To this end, we develop a new framework for solving the three-particle 
evolution equations, dubbed Hamiltonian approach in what follows. A short account of 
the same technique is presented in our letter |]12|, where it was used to calculate l/N^ 
correction to the evolution in flavor-nonsiglet sector. 

The basic idea of our approach can be explained as follows. As it follows from ( |3.8| ), 
the mixing matrices Vat is ( p.4| ) do not have any obvious symmetry and, in general, are 
quite complicated. In particular, they are not hermitian and their eigenvectors are not 
orthogonal to each other. On the other hand, by a numerical diagonalization, Eq. (|3.6| ), 
one flnds that all eigenvalues of these matrices (anomalous dimensions) are real for ar- 
bitrary A^. This property is not obvious and allows to suspect some hidden symmetry 
of the problem, which is not manifest in the particular representation of the evolution 
equations ( |3.4| ) involving only forward matrix elements of the operators. We will argue 
that this symmetry indeed exists and is nothing else as the familiar conformal symmetry 
of the QCD Lagrangian. 

The conformal symmetry manifest itself through the SL{2, M) invar iance of the renor- 
malization group equations describing the evolution of the local twist-3 operators includ- 
ing operators with the total derivatives. The SL{2,'R) symmetry of these evolution 
equations is obscured by the restriction to forward matrix elements of the operators in 
( p.4| ) (or, equivalently, the condition that momentum fractions of the partons sum to zero 
in (|3.1|), cf. ( |2.16|) , ^1 + ^^2 + ^3 = 0). Since the operators containing total derivatives 
have vanishing forward matrix elements, it seems natural to neglect their mixing with 
the twist-3 operators ( p.24| ) in the discussion of deep inelastic scattering. But it is this 



reduction that comphcates the structure of the evolution equations if it is imposed from 
the beginning. 

Our approach relies on the conformal symmetry of the evolution equations and can be 
illustrated by the following scheme indicating a chain of transformations on the matrix 
of the evolution kernels Vn- 

forward \ / non-forward \ / hermitian kernels \ 

non-hermitian kernels J \ hermitian kernels J \m conformal basis J 



Instead of dealing with the non-hermitian "forward" mixing matrices Vn of dimension lq+ 
ig = (A^-|- 1) -|- [A^/2], with iq and ig being the total number of quark-antiquark-gluon and 
three-gluon forward matrix elements, we choose to consider much bigger matrices (but 
hermitian with respect to the so-called conformal scalar product) of dimension [iq{iq — 
1) + ig{ig — l)]/2 that take into account the mixing with the operators containing total 
derivatives. Diagonalizing the thus defined "non-forward" evolution kernels we expand its 
eigenstates over the basis of "spherical harmonics" of the conformal SL{2,M.) group and 
obtain much simpler matrix equation of the coefficients in this expansion (see Appendix 
C). Thanks to the conformal invariance, the corresponding matrix, defining the non- 
forward evolution kernels in the conformal basis representation, has smaller dimension, 
iq + ig, and is now hermitian. We then make a forward projection at the very end. 

3.1 Coefficient functions of local operators 

An arbitrary local three-particle operator O^^a is defined by the set of coefficients in the 
expansion over the standard basis of operators built from the elementary fields $fc and 
covariant derivatives (cf. ( ^.5| )), schematically 

ON,a= Yl ^feifefc3(^'^*l)(^'^$2)(^'^$3), (3.11) 

ki+k2+k3=N 

or, equivalently, by a characteristic homogenous polynomial of three variables 

^nA^,,X2,xs)= Y1 <lkA4'4\ (3.12) 

kx+k2+ks=N 

which we shall refer to as the coefficient function of the operator (not to be mixed with 
the coefficient functions in the OPE). The rationale for the name is that with the help 
of the coefficient function the local operator can be "projected out" of the corresponding 
nonlocal operator, in our case 



O 



zi=0 



(3.13) 



Zi=0 



where \E'^(xj) and \E'^_^(xj) are homogenous polynomials of degree Nggq = N and Nggg = 
N — 1, respectively.^ The construction of the multiplicatively renormahzable operators 

^Tlic difference in the number of derivatives is compensated by the different dimensions of the quark 
and giuon field. 



-I o 



( p.l3| ) is equivalent to finding of the appropriate coefficient functions. The normahzation 
color factors have been included in ( |3.13|) for later convenience (see Eq. (|3.65|) below). 

To expose the symmetries of the problem, it is convenient to start with the evolution 
equations for the corresponding nonlocal operators Sfj,{zi) and 0^{zi) defined in (|2.11| ) 
and (|2.12| ), respectively. They have the following general form 



5 ^2, Zs) 



dQ 

Q'^^7^0^{zi,Z2,Z3] 



dQ 



4:71 

Ait 



Hqq (Zi, Z2, Z^) + Hgg O 

fjL yZi^ Z2i z^j 

Hgq {Zi^ Zi-, Z'i) + Hgg Ofj, {Zi^ Z2, Zs] 



(3.14) 



where Zi stand for the light-cone coordinates of the quarks and gluons, and Hab {a, b = 
q, g) are integral operators acting on quark (gluon) coordinates and describing the inter- 
action between quarks and gluons on the light-cone, HqqS^ (zj) = / dz'j^Hqq{zi\z'j^)S^{z'f^) 
and similar for the other kernels. Note that the short- distance expansion of the nonlocal 
operators and gives rise to the local twist-3 operators ( p.24|) as well as operators 
with the total derivatives. In difference to the previous discussion we do not assume the 
translation invariance, (^i, Z2, z^) ^ {zi + 5,Z2 + 5, ,23 + 5), so that the mixing with 
operators containing total derivatives is included in ( |3.14| ). The explicit expressions for 
the kernels Hab can be found in [^]. They will not be needed for what follows. 

The evolution equation ( |3.14|) has the form of the Schrodinger equation with the 2x2 
matrix of the evolution kernels Hab playing the role of the Hamiltonian. Let '^'^^ ^{zi) 
and '^%_i^{zi) be homogenous polynomials in the light-cone coordinates of quarks and 
gluons of degree N and — 1, respectively, satisfying the Schrodinger equation 



Hqq Hqg 
^91 ^99 





(3.15) 



with the same evolution kernels as in (|3.14| ). The nonlocal operators S^{zi) and 0^{zi) 
can be expanded in terms of these functions with certain operator-valued coefficients 



Sjj,{zi) 
0^{zi)^^^ 



El $!^"^?^ I <^nAQ'), (3.16) 



-f9 

■ N-l,a 




where the subscript in the l.h.s. stands for the normalization scale. It follows from 
the evolution equation that the local operators ON,a that appear in this expansion are 
renormalized multiplicatively and their anomalous dimensions are determined by the 
corresponding energy eigenvalues SN,a 



Tn 



£ 



N,a 



(3.17) 



with the subscript a = 0, 1, [3A^/2] enumerating different solutions to ( [3.151 ). The scale 
dependence of Otv,^ takes, therefore, the standard form 



(3.18) 



-1 A 



(VI/ 



So far we have introduced two different sets of polynomials: (\E' 



N 



^'^_i(zi)) and 

i{xi)). The former set defines the expansion (|3.16|) of nonlocal operators 
S^{zi) and 0^{zi) over the complete set of local multiplicatively renormalizable operators 
On, a) while the latter determines the particular form of the multiplicatively renormal- 
izable operators Oat^o, cf. ( p.ll| ), ( |3.12|) . It follows from the definition ( |3.13|) and ( |3.16| ) 
that 



2iV,vl/^ (9, 



-l,a' 



(3.19) 



and in this sense the polynomials \l/ are dual to the coefficient functions ^ . In what 
follows we shall refer to \E'— functions as coefficient functions of a local operator in the dual 



representation. One can determine the functions ^ and ^ from the orthogonality 

and $f 



condition ( p.l9| ) provided the complete set of eigenfunctions 



N-l,a 



of the 



Schrodinger equation ( |3.16| ) is given. This task seems complicated, but in fact is not. 
Using conformal symmetry, we will be able to find an explicit expression connecting the 
"direct" \E'(xj) and dual "^{zi) coefficient functions of a multiplicatively renormalizable 
operator. The answer is given in Eq. (|3.37| ) below. 



3.2 Conformal symmetry 



A remarkable property of the evolution kernels in (|3.14|) is that they are invariant under 



the projective transformations of the light-cone coordinates of quark and gluons, z^, 

Zk — ^ , ad — be = 1 . (3.20) 

czk + d 

This invariance has its roots in the conformal symmetry of the QCD Lagrangian and the 
transformations (|3.2CI|) form the SL{2, R) (collinear) subgroup of the full conformal group 
acting on the fields "living" on the light-cone. As well known, the conformal symmetry of 
QCD is broken by quantum corrections. However, since the leading-order renormalization 
group equations are driven by tree-level counterterms, they have to respect the symmetry 
of the QCD Lagrangian. 

The action of the SL{2,'R) transformations (|3.2(]| ) on (quantum) fields ^a{z), where 
a = q ,q , g ,g corresponds to q{z), q{z), Gx±{z), Gx±{z), respectively, is defined as 



<l>,{z)-.{cz + d)-'^<^J^^] (3.21) 

\cz + a J 



and is described by three generators L+, L_ and Lq that can be realized as first-order 
differential operators acting on the field coordinates: 

L-$,(Z) = -d^^aiz) , 

L+<l>,(^) = {z^d, + 2jaz) $„(z) , (3.22) 

Ll<l>,{z) = {Zd,+Ja)<^a{z). 

Here ja = {la + Sa)/2 is the conformal spin of the field ^a{z), with being a canonical 
dimension (3/2 for quarks and 2 for gluons) and Sk the spin projection on the light-cone 



direction, Sp^^a 
values Sf, = s 



isa^a- III the case at hand the spin projections have their maximum 



1/2, 



^9 



1, leading to 



J5 = 3/2. 



(3.23) 



In order to unify the notation, we introduce two 2x2 matrices of the evolution kernels Ti. 
and the generators of conformal transformations Ck {k = ±, 0) in the quark- ant iquark- 
gluon and three-gluon channels 



H 



Hqq Hqg 

^gq ^gg 



Tk 

igq 





fk 

^9~99 



(3.24) 



where L^g^ and {k = ±, 0) are the total three-particle SL(2, M) generators acting on 



aaa 



antiquark-quark-gluon and three-gluon coordinates, respectively: 



991? 



Tk 

^999 



The conformal invariance of the evolution equation is stated as 



[n,Ck] = [n,£'] = [c'\c 



■r2 







where 



2^ + 



(3.25) 

(3.26) 
(3.27) 



is the three-particle quadratic Casimir operator. 

Thanks to the conformal invariance ( |3.26 ) the solutions to ( p.l5 ) can be classified 
according to the representations of the SL{2, M) group. Namely, we can impose additional 
constraints on the eigenf unctions 



Cn-J 



0, 

(3.28) 
Since 

the eigenf unctions are homogenous polynomials in the light-cone coordinates Zk, the first 
equation in (|3.28|) is automatically satisfied provided 



where the S'L(2,M) generators are given by the same expressions as in (|3.22 



J 



7 



9 



999 



999 ■ 



(3.29) 



where Nggg and Nggg count the total number of covariant derivatives in the corresponding 
local operator, Nggg = Nggg + 1 = N. Notice that for J = 7/2, or equivalently Nggg = 
there exists no three-gluon contribution. The second condition ensures that \E'^ and 
are invariant under translations of the light-cone coordinates. This requirement defines 
the so-called highest weight of the discrete series representation of the SL{2,M.) group, 
labeled by the conformal spin J = N + 7/2. An infinite tower of solutions to ( ^.151 ) can 
be obtained from the highest weight by the repeated application of the 'step-up' operator 

' (3-30) 



Since £+ commutes with the Hamiltonian, Eq. ( |3.26| ), all states "\I/Ar, with different 
n = 0, 1,2, . . . have the same energy. As we will show in a moment (see Eq. ( p.36| )), 
the corresponding operators ^On are just those obtained from the highest weight state 
operator by adding the n-th power of the total derivative and, therefore, they do not 
survive upon taking a forward matrix element. Note that the expansion in (|3.16|) formally 
includes operators with arbitrary powers of total derivatives, but we can ignore their 



contribution and concentrate on studying the properties of the highest weights (|3.28| ) 
only. 

Going over from the dual coefficient functions "^{zi) to the coefficient functions "^{xi) 
defined in (|3.13| ) corresponds to going over to a different (non-standard) representation 
of the conformal group. Using the relation ( |3.19| ) and requiring 



vi/(5,„a,„a,3) te'°$(^i,^2,^3) 



(3.31) 



one finds the following representation of the SL{2, 
coefficient functions '^{xi,X2,Xs) 



generators on the space of the 



Ll<il{xk) = {xkd^, +Jk)'^ixk), 

L+^{xk) = -a;fe^(xfc), 

Lk ^(a^fc) = {xkdl^ + 2jkd,J^{xk) 



(3.32) 



with the conformal spins jk defined in (|3.23| ). Note that these expressions are more 
complicated compared to the standard expressions ( p.22| ). Eqs. ( |3.32| ) and ( |3.22| ) define 
the SL{2) generators in two different representations and, as such, they are related to 
each other through a transformation 



^{Zi 



[T^izi 



(3.33) 



that maps into each other the coefficient functions in two different representations. The 
explicit form of the T— transformation is given by 



^{z.:) = ^{d,;)\{{i-x.,z.:)-^'^ 



i=l 



n 



r(2j-,) 



e *'^(zA). 



(3.34) 



To verify this relation, use (|3]2|) and (|02|) to check that L^^{zi) = -[T(L^^)](zi) and 
L^'^lzi) = [T{L^'^)]{zi). The conformal constraints ( p.28|) on the coefficient functions 
corresponding to the highest weight look exactly as before: 



[Co - J] 



Xi 



- J( J - 1 



0, (3.35) 



with the SL{2) generators defined in the representation ( |3.32| ). Note that the "step- 
up" operator L4. has become very simple and its action consists of adding the sum of 
derivatives acting on each of the three fields, as we anticipated: 



N 



Xi 



[Xi +X2+ Xg) 



(3.36) 



Finally, applying the transformation ( p. 341 ) to the coefficient function ( |3.12| ) we obtain 
the following expression for the dual coefficient function 



ki+k2+ks=N 



"Jkik2k3 ^1 ^2 ^3 



2ji) r{k2 + 2j2) r(A;3 + 2j3 



r(2ji 



r(2j2 



r(2j3 



(3.37) 



This relation establishes the one-to-one correspondence between the coefficient functions 
( p.l2| ) and their dual counter-parts (|3.37|) . 

Much of the following discussion is based on the fact that the coefficient functions 
of multiplicatively renormalizable operators '^N,a{xi) satisfying the highest weight con- 
dition (|3.35|) are orthogonal with respect to the so-called conformal scalar product. This 
property becomes crucial in establishing the hermiticity of the evolution Hamiltonian in 
( p.l5| ). The hermiticity property will be quite helpful in the further analysis and as we 
argue below is a direct consequence of the conformal symmetry. 

Eq. ( |3.19[ ) suggests to define the following scalar product on the space of coefficient 
functions given by ( p.l2| ) 



(3.38) 



For the coefficient functions of local operators without total derivatives that satisfy the 
constraints ( |3.28| ) (i.e. those that we are interested in) one can equivalently rewrite the 
definition in ( |3.38| ) in a more familiar integral form: 



{'^N,a\'^N,p) 



r(2ji + 2j2 + 2j3 + 2N) 

r(2ji)r(2j2)r(2j3) 



[dx] Xi 



2ii-l 2j2-l 2i3-l 



X 



ix^^i 



(3.39) 

where ju are the conformal spins of the operators entering ( p.ll| ). Here, the integration 
goes over the region < < 1, Xi + 0:2 + 0:3 = 1 and the integration measure is defined 
as 

[dx] = dxidx2dx3 6{xi + X2 + a;3 — 1) . (3.40) 

The SL{2) generators ( |3.32|) are (anti)self-adjoint operators on the space of the coefficient 
functions endowed with the scalar product (|3.38|) . Indeed using Eqs. ( p.31| ), (|3.32|) and 
( p.33| ) it is straightforward to verify that 



As a consequence, the two-particle Casimir operators L^^ 
are the self-adjoint operators. 



{^N,a\Ll^N,/3) = {Ll^N,a\^N,/3) ■ (3.41) 

2 - {Li + Lif defined as 



3.3 Helicity basis 

As we will see in the next section, the evolution Hamiltonians (|3.24|) can be written in 
terms of the two-particle Casimir operators L'jf^ = {Li + in the dual representa- 
tion. This property ensures that the Hamiltonians inherit hermiticity properties of the 
generators and are self-adjoint operators as well. As a consequence, their eigenvalues 



alias the anomalous dimensions (|3.17|) are real and the corresponding eigenfunctions are 
orthogonal to each other 

(3.42) 



In fact, for three-gluon operators there is a complication that the solution to the Schro- 
dinger equation (|3.15| ) has to be found on the subspace of functions \E'^^(zj) that are 
antisymmetric under the exchange of the first and third gluon zi ^ z^. The permutation 
operator P13 that projects onto the states with correct symmetry is not a self-adjoint 
operator so that one has to be careful^. In physical terms, the problem arises because, 
as we will see in a minute, the three-gluon operator contains a sum of contributions 
of gluons with opposite helicity, antisymmetrized because of the crossing symmetry. The 



way out pO[ is therefore to write down the evolution equation for the helicity eigenstates 
and restore the crossing symmetry at the end. 

The construction of the helicity basis is based on the decomposition of the quark and 
gluon fields entering the definition of the nonlocal operators S^{zi) and 0^(2j) into the 
components of different chirality 

qH^) = ^^?(^) ' = ^"(^)^^ ' (^) = I [^M." ± ' (3-43) 

where e^^ = ef^uaisPanfs/ipn). The fields defined in this way satisfy the conditions 

^,q^{z) = ±q^{z) , ^ei,Gt{z) = ±G^{z) (3.44) 

and describe quark, antiquark and gluon of a definite helicity. Making this decomposition 
one obtains the expressions for the nonlocal operators Sf^{zi) and O^(zj) in terms of the 
chiral fields. We remind that /i = 1, 2 is a transverse index and in order to construct 
the three-particle states with definite overall helicity one has to take particular linear 
combinations, projecting onto the two complex vectors in the transverse plane: 

= + '^±1 ' = - ■ (3.45) 

We find for the quark-antiquark-gluon operator 

S^{zi,Z2,Zs) = {zi)G:l^{z2) fiq^ {Z3) + q'^ {z3)G^{z2) flq' (zi) , 

Sj;j{zi, Z2, Z3) = q+{zi)G^{z2) fiq'izs) + q'{z3)G^{z2) fiq'^izi) . (3.46) 

Notice that the quark and the antiquark have opposite helicity in both cases, and helicity 
of the gluon ±1 concides with the total helicity of the system. (We tacitly imply that the 
momenta of the three partons are aligned along the same light-cone direction defined by 
the proton momemtum p.) The similar decomposition of the three-gluon operator looks 
as follows 

Owizi, Z2, Z3) = - [Tu,{zi, Z2, Z3) — Tyj^z^, Z2, Zi)] , 

Otu{Zi, Z2, Z3) = ^ [Tj;j{Zi, Z2, Z3) - Tj;j{z3, Z2, Zi)] , (3.47) 



where the notation was introduced 



Tw{zi, Z2, Z3) — -^f"'^'^G% (2;i)G'^^(z2)G'^^(z3) , 
Tw(^i, Z2. Z3) = '-lr'^Gl-^{z,)G'^~{z2)G%-{z3) . (3.48) 



^The situation is in fact very similar to the Hartree-Fock construction of the completely antisymmetric 
fcrmionic wave functions in quantum mechanics. 



The operators Tw{zi) and T^^ji^Zi) describe the state of three gluons with the total hehcity 
±1, respectively. 

It follows from ( p.47|) and ( p.48|) that the operators 0^(zi, ^2, -23) and T^(zi, 22, ^3) are 
antisymmetric under the interchange of gluons with the same and the opposite helicity, 
respectively: 

0^(^1,^2,2:3) = —0^(^3,^2,2:1) = —P3iOw{zi^Z2-,Z3)^ 

T^{zi, Z2, z-i) = -T^{zi, Z3, Z2) = -P23T^{zi, Z2, Z3) (3.49) 
with Pik being the permutation operators. Using this, we can invert (p.47|) to get 

Tw{zi, Z2, Z3) = Ow{Zi, Z2, Z3) + Ow{z2, Z3, Zi) — Ow{z3, Zi, Z2) 

= (1 + P12P31 - P23P31) 0^(^i, ^2, Z3) . (3.50) 

The same relation holds between the operators T^j and Ow 

The relations ( p.47| ) and (|3.50|) allow one to rewrite the evolution equation ( p.l4| ) for 
the three-gluon operator 0{zi) in terms of T{zi) with definite helicity. The difference 
amounts to the following redefinition of the evolution kernels ( p. 24 ) 7i — »• 
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Hqh 




\ Hhq 


Hhh J 
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0(1 + Pi2 - P23 




Hqq Hqg 

Hgq Hgg 




Oi(l-P 



31 



(3.51) 



Note that the kernel Hqq is not affected by this transformation. Despite the fact that the 
two evolution equations are obviously equivalent, the premium in dealing with helicity 
operators Tw{zi) (or Tyj{zi)) is that, as we will show below, the evolution kernel Ti^ 
becomes hermitian on the space of the coefficient functions. The original kernel in (|3.24| ) 
is not hermitian due to the presence of additional permutation operators in ( |3.51| ). The 
explicit expressions for the operator Ti^ will be given below. 

According to ( |3.47| ), going over from 0{zi) to the helicity states T{zi) is equivalent to 
the following ansatz for the dual coefficient function of the three-gluon operator: 



-^%_^{Zi,Z2,Z3) = - 



^N-l(^l, Z2, Z3) - ^%-i{z3, Z2, Zi) 



(3.52) 



where the coefficient functions '^%_i{zi) are defined by the short distance expansion of 
T^(2;j). Applying the equivalence relation (|3.50|) one finds 

^%_iizi, Z2, Z3) = %_i(2;i, Z2, Z3) + ^%_^{z2, Z3, zi) - %_i(2;3, zi, Z2) . (3.53) 



Substituting ( p.47| ) into ( |3.13|) , one rewrites the multiplicatively renormalizable oper- 
ators as 

0^,„= [2iV,v^^ J9,J5^(;^,) + n^vI/^_i_J9,jT^(z,)]|^^^^ , (3.54) 
and obtains similar relations between thus defined coefficient functions '^^_i^{xi) and 



We find that the local three-gluon operators are projected out of the T{zi 



by the following coefficient function 



*jv-i,a(a^i> 2:2, 2:3) = - [^^_i^„(a;i, X2, X3) - ^^^.^ „(xi, X3, Xs)] 



(3.55) 



that has the same symmetry ( |3.49| ) as the hehcity operator itself: 



^5^r-i(a;i,a;3,X2) = P23 X2, 0:3) = -^5^_i(a;i, Xs, X3) . 

The inverse relation looks as 



(3.56) 



N-l,a 



(xi, X2, X3) = ^(xi, X2, X3) + ^?^_i_„(x3, xi, X2) - ^5^^-1,0(3^2, 3:3, Xi) . (3.57) 



Notice that the relations between the two sets of the coefficient functions are different in 
the direct and dual representations. Nevertheless, it follows from ( p. 521) and ( p.57|) that 
the coefficient functions '$%_^{zi) and \l/5^r_i(xj) satisfy the same conformal constraints 
( p.28| ) and (|3.35|) as \l/^_i(zi) and \l/^_]^(xi), respectively. 

Last but not least, the original Schrodinger equation ( p. 151) can be transformed into 
a Schrodinger equation for the coefficient functions \E'^(xi) and \E'^_^(xj) 



Hqq Hqh 
Hhq Hhh 



N 



(3.58) 



where the evolution kernels "without a hat" in the helicity basis, if^^, are obtained from 
the corresponding kernels if^^, ( |3.51| ) through the T— transformation ( pi 



Hab = T„ ^Hnt^. 



ab^b 



(3.59) 



3.4 Conformal basis 

Solution of the evolution equations (|3.58|) can be simplified significantly by expanding 
the coefficient functions \l/^(xj) and \l/^_^(xj) over the basis of "spherical harmonics" 
consistent with conformal symmetry constraints ( |3.28| ). 



3.4.1 General construction 



For a generic three-particle operator, (|3.11| ) and ( p.l2|) , such a 'conformal basis' can be 
constructed as follows |2T|. Conformal symmetry allows to fix the total three-particle 
conformal spin J = ji + J2 + is + ^ of the state that translates to the condition that 



the corresponding coefficient function satisfies Eqs ( |3.28| ). We define the set of functions 



yj^^^^ by requiring that they obey Eq. (|3.28|) and, in addition, have a definite value of 



J 3 

the conformal spin in one of the two-particle channels, (31) for definiteness. The latter 
condition leads to 

Ll,Yfl^\x,) = j{j - l)Yif\x,) , (3.60) 
where j = ji + js + n and n = 0, A^. 



Taken together, Eqs. (|3.60|) , ( p.28| ) determine the functions Yj'^-^^'^ uniquely and have 



the following solution: 



) = rjj{xi +X2 + xsY'^'^^xi + x-sY'^''^' 



p(2j2-l,2j-l) 



Xi — X2 -|- X3 ^ p(2j3-l,2ji-l) I Xi ~ X3 



X1 + X2 + X3 



Xi +Xs 



(3.61) 



r» -1 



Here P^'^\x) is the Jacobi polynomial and rjj an arbitrary normalization factor. The 



basis functions Yj'f^^'^{ orthogonal with respect to the scalar product ( |3.39D . Re- 



quiring that = 1 we fix the normalization to be 



r 



-2 



r(2J) r(j+ji-j3)r(j-ji+j3) 



r(2ji)r(2j2)r(2j3) r(j - ji - ja + i)r(j + ji + ja - i)(2j - 1) 

r(j-j+j2)r(j + j-j2) 
r( J - J - j2 + i)r( J + J + j2 - 1)(2 J - 1) ■ 



(3.62) 



The above construction of the conformal basis involves an obvious ambiguity in which 
order the spins of partons are coupled to the total spin J. Choosing in ( p.60| ) a different 
two-particle channel one obtains a different conformal basis related to the original one 
through the matrix Vt of Racah 6j-symbols of the S'L(2,R) group 

yjT\^i)= E ^AJ)y!,f\^i)- (3-63) 

Properties of the Racah 6j-symbols as well as their explicit expressions in terms of the 
generalized hypergeometric series 4-^3(1) are summarized in Appendix A. 

3.4.2 Quark-gluon conformal basis 

We now specify the above general construction to the particular cases of quark-antiquark- 
gluon and three-gluon operators and define 

^7V,fc(^l5^2,a;3) = ^7V+7/2.A:+2(^l'^2,a;3) . 

J1=J3=1 j2=3/2 

1^7v_i,fc_i(xi,X2,X3) = yi+J/2,n+2(a^i'2;2,a;3) , (3.64) 

Jl=j2=j3=3/2 

where the prefactors have been chosen for later convenience. The three-particle conformal 
spin is given by J = + 7/2 = jq + jg + jq + N = jg + jg + jg + N — 1 in both cases, 
and the conformal spin in the (31)— subchannel is equal to j = k + 2 = jq+jq + k = 
jg + jg + k — 1. Notice that j > 3 and, therefore, the basis functions Yj^ f^_^{xi) in the 
gluon-gluon subchannel are well-defined for k — 1 > 0. In what follows we assume that 

Each of the two sets of functions in ( |3.64| ) forms an orthonormal basis with respect to 
the conformal scalar product ( |3.39|) whose explicit form depends on the conformal spin 



of the particles and, therefore, is different for the quark-antiquark-gluon and three-gluon 
systems. This suggests to define the scalar product on the space of two-dimensional 
vectors of coefficient functions in the following form 

(*i|^2) = 2N,{n'^l) + n;(vl/t|vl/^) , vl>, = , (3.65) 



where the scalar product in each sector, (\E'^|\E'2) and (\I/^|\E'2), is obtained from (|3.39| ) 



by substituting ji = J3 = 1, j2 = 3/2 and ji = j2 = ja = 3/2, respectively. As we show 
below such choice of the scalar product ensures the hermiticity of the evolution kernels. 



We shall look for the solutions to the evolution equation ( |3.58| ) in the following form 

\ 

(3.66) 



TV 



AT / "^>fc yi (j.. 



E 

fc=0 



"iv.fc yh I 



J 



with and 



u 



N,k 



being the expansion coefficients. In this way, the evolution kernels 
finally become symmetric and real matrices acting on the vector of the expansion coeffi- 
cients. We fix the normalization of the coefficients by requiring the eigenstates (|3.66|) to 
have the unit norm 



TV 



where u%q 



11*1 
0. 



u 



7V,fcl 



+ U 



N,k\ 



(3.67) 



k=0 



3.5 QCD evolution kernels 



To one-loop accuracy, the evolution kernels are given pOf by the sum of two-particle 
Hamiltonians describing the pair-wise interaction between quarks and gluons on the light- 
cone: H = Hi2 + H23 + H31, or, equivalently H = H12 + H23 + H31. Conformal invariance 
( p.26| ) then implies that each pair-wise Hamiltonian Hik only depends on the sum of 
conformal spins of the interacting particles, e.g. in the "direct" representation Hik = 
H{Jik), where 

= {U + Lkf = MJ^k - 1) (3.68) 

and the SL{2) generators (a = ±, 0) are given in Eq. ( p. 221 ). The explicit form of this 
dependence can most easily be obtained by comparing the eigenvalues. To this end it 
is sufficient to calculate the one-gluon exchange diagrams in a simplified situation when 
momenta of the contributing partons sum up to zero. Alternatively, one can start with 
the known integral representation for the evolution Hamiltonian ^ and project it onto 
the conformal basis ( |3.64 ). 



3.5.1 Diagonal evolution kernels 

The diagonal quark evolution kernel is given by 



-^qgl^ N,n/ 



2ni 



l-^ N,n/ 



(3.69) 



where b = ll/3Nc + 2/3n/ is the lowest-order coefficient of the QCD /3— function. The 
first term in brackets stands for the flavour-nonsinglet evolution kernel (see below) and 
the second term proportional to 5ji3,2 comes from the additional Feynman diagram in 
which the quark and the antiquark annihilate to produce a gluon that splits again into 
the gg-pair. Since the quark and the antiquark are produced in this way in one spatial 
point, contribution of this diagram is different from zero only for J13 = jg + jg = 2. The 
flavour-nonsinglet Hamiltonian if 5+ can be represented as pOl E2|, E3, E4 



Hs+ =NrH^''^ 



(3.70) 



where 

H^'^ = V^fiJi2) + U!^^\j,s) (3.71) 

H^'^ = V,^i\Ji2) + U^l\j,s) + C/W(Ji3) . (3.72) 

Here, the notation was introduced for the two-particle quark-quark and quark-gluon 
kernels 

V^f (J) = ^{J + 3/2) + V^( J - 3/2) - 2^{1) - 3/4 , (3.73) 
[/(;)( J) = V'( J + 1/2) + ^l^iJ - 1/2) - 2V'(1) - 3/4 , 

^a)(n. (-^y-'^' C;a)(n._izit!^l 

(j-3/2)(J-l/2)(J+l/2) ' 2(J-l/2)' 

t^i?(^) = ^[V'(^-l) + V'(J+l)]-V'(l)-3/4, 

where xp{x) = dlnT{x)/dx. 

The diagonal gluon kernel is defined as 



where 



Hhh\Y}^-,,k-i) = [Hz/2 - V^,2] (3.74) 

= 2 [[/,,( J12) + t/,,( J23) + f/..( J31)] - W , (3.75) 
V„2^V,,{h,)+V,M (3.76) 



and 



U,,{3)^m-m, (3.77) 
2 3(1 + (-1)^-) 

j(j-l) (j-2)(j-l)j(j + l) 



V..(J) = 777^ + ;\:7^^^^ .^ ■ (3-78) 



We notice that the kernel 11^^ is invariant under permutations of the two gluons of the 
same helicity 

[Hhh,P2z]^^. (3.79) 

Note that the pair-wise diagonal Hamiltonians Hik and Hi^ have the same functional 
dependence on the Casimir operators, i.e. — h{Jik) and Hih — h{Jik) with the same 
function h. 

3.5.2 Off-diagonal evolution kernels 

The off-diagonal kernels describe the mixing between quark-antiquark-gluon and three- 
gluon states. The conformal symmetry implies that the conformal spin of both states 
should be the same, and this appplies both to the total conformal spin of the three-parton 
system and the conformal spin of the parton pair involved in the mixing. It follows that, 
therefore, acting on the basis function in the quark sector, \Y^i^), the evolution kernel 
Hhq transforms it into the basis function in the gluon sector, \Yl^_i^f,_i), with the same 



A^, k. Similarly, the evolution kernel Hqh niaps gluon conformal basis into the quark one. 
As a consequence, the off-diagonal kernels can be written down in the following form 



H, 



qh 



23 



31 > 



l-B. 



23 



hq 



31J 



(3.80) 



where the operators Wqh and Whq is defined as 



W, 



qh 



\^ N- 



l,k- 



-l) 



7^3- Ji3 + 2(-l) 



Jl3 



MJl3 - l)v/(^13 + l)(^13-2) 



'^13 



Ji3 + 2(-l)-^« 



Jl3(^13 - l)v/(^13 + l)(^13-2) 



\Y, 



h 



(3.81) 



The following comments are in order. The off-diagonal kernels Hgh and H^q originate 
from the Feynman diagrams in which quark and antiquark annihilate into two gluons 
of opposite helicity. As a consequence, these kernels depend on the conformal spin in 
the quark-antiquark subchannel, J13, and the additional projector (1 — P23)/2 takes into 
account antisymmetry of the three-gluon state under the interchange of gluons with the 
same helicity, Eq. p.49| ). 

It follows from ( |3.81| ) that the off-diagonal kernels are adjoint to each other, Hqg = 
H^q, with respect to the scalar product ( |3.65| ). Indeed, calculating the matrix elements 



of the off-diagonal kernels between the states and defined in ( |3.66|) and taking 

into account their symmetry properties (|3.56|) one finds 



h \ 
N~l/ 

N 



(fc + l)(fc + 2) + 2(-l)^ 



K,ky< 



Af-l,fc-l 5 



(3.82) 



where we have substituted J31 = 
elements of the diagonal kernels. 



Hi 



Hqg and Hl^^ 



k + 2. Repeating the similar calculation for the matrix 
and {'^'^_-^^\Hhh'^%-i) , one makes sure that 
Hhh- We conclude that the matrix of the evolution kernels entering 



the Schrodinger equation (|3.58|) is a hermitian operator on the space of the coefficient 
functions endowed with the scalar product ( ^.65| ). As a consequence, its eigenvalues £n 
are real and the corresponding eigenfunctions are orthogonal to each other. 



3.6 Expansion in conformal operators 

Once the Schrodinger equation in the helicity basis ( p.58| ) is solved, we can easily restore 
the gluon coefficient function ^ using the symmetry relation (|3.57|) and reconstruct the 
multiplicatively renormalizable operators ( p.l3| ). Their reduced forward matrix elements 
can be expressed in terms of the multiparton distributions as follows: 



1^ 



N,a\ 



(3.83) 

Notice that in contrast with ( |3.39|) the integration goes here over the region ^1+^2+^3 = 0. 
Comparing this expression with the expansion is the standard operator basis ( p.5|) and 



assuming that the eigenvectors '^N,a is normahzed to unity one gets the equivalence 
relations 



N~k 

0<k<N 



nj K-lJ^^) „ , = E ^-(^) i^T^t'""" - ef-'-'"C) • (3.84) 

^ 0<m<[JV/2]-l 

Next, expanding the OPE coefficient function $^(^1,^3) over the eigenstates of 
the evolution kernels 



J . v%^'^-i,«te 



(3.85) 



we can decompose moments of the structure function g2{x,Q'^) ( ^.31 ) in multiplicatively 
renormalizable contributions ( p.83| ) as 



4 



^72'L(iV + 3,g2) = -(e;)^ 5^ w^AOnAQ')))- (3.86) 



a=0 



The expansion coefficients WN,a can be expressed in terms of C^k^N) and C^^{N) by 
comparing the coefficients in front of different powers of .^1 and .^3 in the both sides of 

(S). 



3.6.1 Quark coefRcient function 

There exists, however, a more efficient way of finding the same coefficients. To this end, 
observe that the OPE coefficient function $^^3(^1, —^1—^3, ^3) can uniquely be continued 
from the hyperplane ^1 + ^2 + ^3 = to arbitrary values of by requiring that the thus 
defined function $^_,_3(^i, ^2, ^3) satisfies the conformal constraints (|3.35| ) (and coincides 



with $^_,_3(^i, —,^1—^3, ^3) at ^2 = — 'Ci — 'Cs). To find an explicit expression, let us expand 
^Af+3(^i' ^2, ^3) over the conformal basis ( |3.64| ) 



N 



$^+3(xi, X2, X3) = E (t)%^k K,ki^i) (3.87) 



k=0 



with = {^N+^iy^ k) ■ Projecting out the hyperplane ^jXj = 0, we find using ( |3.61| ) 
that the basis functions are reduced to 



^N,ki''^ij 



~ (x, + x.fPt'''^ ( ) (3.^ 



and form an orthogonal basis on the subspace Xi + X3 = 1 and < Xi,X3 < 1. This 
property allows us to calculate the expansion coefficients 07v,fc as 

4>Nk^ dxidx3 6{l - xi - X3) xiXsP^^'^H ^^ '^^+3(3^1, -xi - X3,X3) . (3.89) 

Jo \Xl + Xs/ 



Substituting the actual expression for ^^^^{^i) ( |2.22| ) and performing the integration one 
arrives at 



^ ( {k + l){k + 2){2k + 3){N-k + l){N-k + 2) T^N + 3) ^ _ 



8{N + 3){N + k + A){N + k + 5) r(2A^ + 6) 

Using these coefficients one can easily calculate the norm of the quark coefficient function 

N 



|<J)5 ||2 



k=0 



'''N,k) 



l + -(l-4V'H-47£;)-- 



(3.91) 



where in the r.h.s. n = N + 3. 

Finally, using ( p. 87] ) and ( p.90|) one finds the expansion coefficients for moments of 
the structure function 



1 ^ 

^ ^ k=0 



9 

N,k 



(3.92) 



The coefficients m^^^ stand for the expansion of a multiplicatively renormalizable operator 
in the conformal basis, cf. Eqs. (|3.66| ) and (|3.67|) . Their dependence on the particular 
operator (index a) is tacitly assumed. 



3.6.2 Gluon coefRcient functions 

The similar procedure can be used to calculate the expansion coefficients of the gluon 
coefficient functions $^ and defined in ( |2.19| ) and ( |2.23| ). The corresponding conformal 
harmonics in the helicity basis are given in the forward direction by 



^N-l,ki^i 



(3.93) 



In order to decompose the gluon coefficient function ( |2.23| ) over this basis we first con- 
struct the same coefficient function in the helicity representation 



1 



(3.94) 



so that P23^%^.;^{xi) 



and then define the expansion coefficients ((>%_i ^ as 



N-l 



-l,k^N -l,k{^i) ■ 



(3.95) 



fe=0 



Using orthogonality of the Jacobi polynomials we can calculate these coefficients as 



dxidx^Sil-Xi-x^) {xix^f P]^''^^ [ "^^ )<l>5^_^3(xi, -Xi-X3,X3) , (3.96) 



r\i-r 



Substituting the gluon coefficient function ( |2.23| ) into ( |3.94| ) and ( p.96| ) and going through 
the calculation one finds the following explicit expressions for the expansion coefficients 
for even conformal spins N 



for even k and 



(A^ + A; + 5)(A^ + 7 + 5/2 A; + 1/2 A;^ 
{k + l){N-k + l) 

{k + 4){N + k + 6) 



h 



N,k 



'^N-l,k 



N,k 



2{N -k) 

for odd n. Here, the normalization constant /iat ^ is given by 



'''N,k 



{k + 1){2 k + 5){N -k){N -k + 1){N + 3)2 r^(A^ + 3) 
8{k + 3){k + 4){k + 2){N + 5 + k){N + k + Q) T{2N + 6) 

Using (|3.95| ) we can calculate the norm of the gluon coefficient function 



(3.97) 



(3.98) 



(3.99) 



N-1 
fc=0 



r^(n) n\n-l] 
r(2n) 16 



■ 2(n-4)(n+ 1) , , , , 



(3.100) 



where n = + 3, as above. 

The second gluon coefficient function, appears in the expression for the mo- 

ments of the structure function ( |2.19| ) due to mixing between quark-antiquark-gluon and 
three-gluon distribution amplitudes. Similar to (|3.94|) and ( p.95| ), one can transform this 
function into the helicity representation and decompose it over the conformal basis 



- [f^^+3(xi, X2, Xs) - fi?^+3(xi, X3, X2)] = ^N-l,kyN-l,kiXi) ■ (3.101) 



A:=0 



Substituting by its explicit expression (|2.23| ) it becomes straightforward to calcu- 

late the coefficients ^ using (|3.96| ). The resulting explicit expression is cumbersome 

and will not be displayed here. Comparing the expansions ( |3.95| ) and ( |3.101| ) one finds, 
however, |]T^ that the two coefficient functions $^(^i) and are to a good accuracy 

proportional to each other 

fi^/(e.)^c(n)$^je.)- (3.102) 

Applying the transformation ( p.94| ) and ( |3.101[ ) to the both sides of this relation and 
using the orthogonality of the conformal basis we can calculate the coefficient c(n) as 



N-l 
ft II 2 ih 



^JV-l,fe^Af-l,A;/ 



1$ 



hi\2 



(3.103) 



k=0 



with the norm ||$^|| given in ( 3.100| ). Going through the calculation one finds that c{n) 
is given at large n by the following expansion: 



c(n) = 1 + — [41nn + 47E-6] + C(l/n3). 



(3.104) 
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Figure 1: The spectrum of anomalous dimensions of flavor-singlet twist-3 operators with 
the mixing between quark-antiquark-gluon (crosses) and three-gluon operators (open circles) 
switched off, cf. (|4l|). 



4 Results 



In this section we present detailed results on the solution of the Schrodinger equation in 
Eq. (|3.58|) . The main advantage of the Hamiltonian approach described in the previous 
section is that it allows to understand qualitative features of the solutions using the 
intuition and a wealth of analytical tools well known from quantum mechanics. For this 
analysis, we decompose the full Hamiltonian in (^.58|) in two parts: 



Tin 



Hg, 



V 



(4.1) 



Hgh 

HhhJ' " yHhg 

The Hamiltonian Tio governs the scale-dependence separately in the 'quark' and 'gluon' 
sectors, whereas the off-diagonal kernel V describes the mixing between the two sectors. 
Our strategy throughout this section will be first to solve the Schrodinger equation for 
Tio and then examine the deformation of the spectrum induced by V that we, somewhat 
imprecisely, refer to as the quark-gluon mixing. The rationale for such a two-step proce- 
dure is that properties of the individual Hamiltonians Hgg and Hhh have been studied in 



detail in recent works |]22|, ^ Q using their relation to completely integrable models. 
We will, therefore, be able to use these results. 



4.1 Diagonal evolution kernels 

The spectrum of eigenvalues of Hq is obviously given by the superposition of the two 
independent spectra of quark-antiquark-gluon and three-gluon operators 



TT .tMo) 
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hh'^ N-1,13 ~ '-'N,I3^ N-l,f3 



and the corresponding eigenfunctions are given by 
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Figure 2: The dependence of the anomalous dimensions of flavor-singlet quark-antiquark-gluon 
(crosses) and three-gluon operators (open circles) on the number of light quark flavors Uf. The 
dotted line corresponds to the linear dependence ~ (2/3)nj, see text. The quark-gluon mixing 



is switched off, cf. (|4.1|) . 



The superscript (0) stands to remind that the off-diagonal mixing terms are omitted 
and the subscripts a and P enumerate different 'quark' and 'gluon' eigenstates corre- 
sponding to multiplicatively renormalizable operators with N (quark) or — 1 (gluon) 
covariant derivatives with the same canonical dimension (and the same conformal spin 
J = N + 7/2). In addition, we require that the 'gluon' eigenfunctions satisfy the sym- 
metry property (|3.56| ). For a given the total number of the 'quark' and 'gluon' eigen- 
states is equal to £q = N + 1 and £g = [A^/2], respectively, so that a = 0, . . . ,N and 
f5 = 0,...,[N/2]-l. 

The results of the numerical calculation of the spectrum for < 20 and nf = 3 are 
shown in Fig. |I[ Note that the gluon eigenstates (open circles) and the quark eigenstates 
(crosses) occupy two bands that lie on the top of each other. For large N the eigenvalues 
(anomalous dimensions) rise logarithmically 



ACplnN < 



< ANJnN, 



4N^ \nN < Eff^p < QN^ In N 



(4.4) 



and the coefficients in front of In are related to the color charges of the corresponding 
(classical) parton configurations. Since ~ 0{N) levels have to fit within the band- width 
~ 0{\n N) ( [4.4|) , the distance between the neighboring levels in general goes to zero. The 



analysis of the 'fine structure' of the spectra p3| , |24| reveals, however, that in the limit 
N ^ (X) three levels remain separated from the rest of the spectrum by a finite gap. 
These three special levels are: the lowest quark level, the highest quark level and the 
lowest gluon level; they will play a decisive role in what follows. 



4.1.1 The highest quark-antiquark-gluon state 

The eigenvalues S'^ ^ and Sf^ ^ depend on the number of light quark flavors Uf. This 
dependence is shown in Fig. ^ for A^ = 8 and < < 10 and reveals the following 
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Figure 3: The coefficients of the expansion ( |4.6D of fc=o(^«) eigenstates of the quark 

Hamiltonian Hqq. N = 8 and = 20 for the left and the right panel, respectively. 



remarkable pattern: In the gluon sector the nj— dependence of all energy levels is linear, 



s ~ 2n//3, and it can be traced to the additive 6— correction to (|3.75| ). In contrast 



to this, in the quark sector all energy levels except the highest one vary very slowly with 
Hf. At the same time, the n/— dependence of the highest quark level is almost identical 
to that of the gluon levels, ~ 2n//3. To understand this property, notice that 

the n/— dependence of the diagonal quark kernel in Eq. (|3.69|) comes entirely from the 
annihilation term (2n//3)5j^3^2- At large Uf this term dominates and the corresponding 
eigenstate is given by 

<i''^^{x,) = Yl,^,{x,) + 0{l/nf). (4.5) 

Remarkably enough, this relation holds true with high accuracy for small values of Uf 
as well, including Uf ^ 0. To illustrate this, we show in Fig. |^ the coefficients in the 
expansion of f^^Q^Xi) over the complete set of the quark eigenstates: 
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a=0 



Ivl/ 



(4.6) 



for 71/ = 3 and two different values of N. Note that the normalization is such that 
ll^/i/^Q II ~ ^ ^ consequence, the sum of the coefficients squared is equal to one. It 

is seen that the overlap with the exact wave function of the highest 'quark' state a = N 
is very close to unity. 

As a yet another test of the approximation for the wave function ([4.5| ) one evaluates 
the corresponding energy (see Appendix A for details) 
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The few first terms of the expansion of the matrix element at large are 
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= 2 


10.933 


10.987 


18.107 


17.993 


17.350 


17.350 




= 4 


12.629 


12.644 


22.587 


22.395 


22.087 


22.160 


N 


= 6 


13.938 


13.946 


25.785 


25.561 


25.474 


25.541 


N 


= 8 


14.994 


14.999 


28.285 


28.046 


28.107 


28.162 


N 


= 10 


15.876 


15.881 


30.343 


30.094 


30.260 


30.304 


N 


= 30 


20.824 


20.828 


41.634 


41.367 


41.847 


41.878 



Table 1: Exact numerical results for the energy of the three 'special' levels (see text) com- 
pared with the calculation using the approximations for the corresponding eigenfunctions. 
The results for the lowest quark-antiquark-gluon, the highest quark-antiquark-gluon and 
the lowest three-gluon eigenstates are shown in the two left, two middle and two right 
columns, respectively. 



The approximation ( [4.7|) is compared with the exact numerical calculations of the energy 
in Table |l| (two middle columns). The accuracy turns out to be very good, better than 
1% for all N. 

Last but not least, we note that according to ( |3.88| ) the wave function entering ( ^4.5| ), 



^Noi^i) ~ (^1 +^3)^ = {^^2)^ , depends on a single (gluon) momentum fraction only 
and, therefore, it defines a local operator ( p.ll ) and ( p.l2|) that contains derivatives 



acting on the gluon but not on the quark fields. All such operators can be obtained from 
a Tailor expansion of the nonlocal operator ( 2.11| ) in which the quark and the antiquark 



are located at the same space-time point and which can be rewritten as a two-gluon 
operator using the equations of motion: 



5^(0, V, 0) = tm9G(,.{vx) MO) = -tG%{vx)DfGl{0) . (4.9) 

Here a, b are color indices and is the covariant derivative in the adjoint representation. 
Thus, to a good accuracy, the quark-antiquark-gluon three-particle operator with the 
highest anomalous dimension is in fact a two-particle gluon operator! 



4.1.2 The lowest quark-antiquark-gluon state 

Since eigenfunctions of the hermitian Hamiltonian Hgg are orthogonal to each other and 
since the eigenfunction corresponding to the highest anomalous dimension turns out to 
be very close to the eigenfunction of the nj-dependent contribution of the annihilation 
diagram that gives rise to the second term in the Hamiltonian ( p.69| ), it follows that the 
wave functions of all other levels have a negligible overlap with this term and, therefore, 
are to a good accuracy nj-independent, in agreement with Fig. ^. 

This means that within the accuracy of ( [4.5|) , the eigenfunctions and eigenvalues of 
all quark levels other than the highest one coincide with those of the fiavor-nonsinglet 



quark Hamiltonian Hs+ that was studied in |]21| , [23| , ^ . In particular, it has been shown 



that the wave function corresponding to the lowest anomalous dimension coincides, in the 
large- A^c limit, with the tree- level coefficient function in the OPE ( |2.22| ) continued from 
the hyperplane ^1 + ^^2 + 'Ca = to arbitrary values of using the conformal symmetry. 
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Figure 4: The coefficients of tlie expansion of $^_|_3(xi)/||$^_|_3|| ( |4.10| ) over the eigenstates 
^T^^a of the quark Hamiltonian Hqq. 
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(4.10) 



The additional factor in the r.h.s. takes into account the different normalization of the 
states. In Fig. ||we show the coefficients of the expansion of '^'Ar4.3(a^j)/||'^*Ar+3ll '^^^^ 
eigenstates of the quark Hamiltonian Hqq for N = 8 and = 20. It is seen that the 
expansion is dominated by the lowest quark state a = 0. The corresponding eigenvalue 
(anomalous dimension) can be calculated as (see Appendix A for details) 



^7V,0 — 



($9 



N+3 



H, 



(0) + 4^ 



N 



Nr. 



(1) 

N 



(4.11) 



where HTD 



and 12 



Ef = 2*(A^ + 3) + 



E 



(1) 

N 



ln(iV + 3)+7K 



4 6 



(4.12) 



(4.13) 



(AT + 3)2 

The quality of this approximation is illustrated in Table |I] (compare the numbers in the 
first two columns). 



4.1.3 The lowest three-gluon state 

Last but not least, we have to work out an approximate description of the lowest three- 
gluon state. As noticed in [Q, the corresponding eigenfunction appears to be close to the 
one-loop gluon coefficient function $^(^i) (|2.23| ) transformed to the helicity representation 
( ^.94| ) (and continued to arbitrary values of C,i using the conformal symmetry): 

*^-!./^=o ^ = ^^^(n+3(^i,^2,X3) -n+3(^i,^3,a:2)) . (4.14) 
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In difference to the above, this approximation cannot be justified in any formal hmit, but 
is no less accurate, as illustrated in Fig. ^ where we show the corresponding (normalized) 
expansion coefficients in the basis of exact three-gluon eigenstates. 

The corresponding approximation for the energy can be calculated as (see Appendix 
A for details) 
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It is compared with the exact result in Table || (the two last columns). Note that the 
energies of the lowest gluon and the highest quark states are very close to each other, as 
seen also from Fig. |l|: 

In fact, the difference is so small that one can treat these two levels as degenerate for 
most purposes. 



4.2 The quark-gluon mixing 

The mixing of quark-antiquark gluon and three-gluon operators is governed by the off- 



diagonal part V of the Hamiltonian ( |4.1| ). It turns out that the mixing is generally rather 
weak and for most practical purposes can be taken into account using the standard 
perturbation theory. To the leading order, the mixing-induced corrections to the 'pure' 



n A 





£n,o 


£n,o — Sn,o 


£-N,N 




£n,n+i 




= 2 


10.739 


-0.193 


17.679 


-0.428 


18.171 


0.821 


= 4 


12.571 


-0.058 


22.324 


-0.262 


22.648 


0.561 


N = 6 


13.912 


-0.027 


25.637 


-0.148 


25.872 


0.397 


N = 8 


14.980 


-0.015 


28.195 


-0.090 


28.420 


0.312 


N =10 


15.867 


-0.009 


30.278 


-0.065 


30.530 


0.270 


= 30 


20.823 


-0.001 


41.600 


-0.034 


42.056 


0.209 



Table 2: The exact energy of the three 'special' levels (see text) with the mixing-induced 
corrections taken into account compared with the energy of the corresponding eigenstates 
with the mixing effects neglected. The results for the lowest quark- ant iquark-gluon, the 
highest quark-antiquark-gluon and the lowest three-gluon eigenstates are shown in the 
two left, two middle and two right columns, respectively. 



quark or gluon eigenstates (^3)' (113) equal to 

E f^^.(o,° X, , (4.18) 



where the notation was introduced for the mixing matrix 

-MO) 

■ N~l,f3 



It is easy to see that (^5v^|V|^S5^) = and, therefore, the energy eigenvalues do not 
receive any corrections to the first order of the perturbative expansion. This explains 
why the mixing-induced corrections to the anomalous dimensions are very small (see 
Table The numerical results presented below are in all cases based on the exact 
diagonalization of the mixing matrix, and the first-order expressions in ( [4.18| ) only serve 
to illustrate the picture. Note that the perturbation theory cannot be used to describe 
the mixing betwen the highest quark-antiquark gluon and the lowest three-gluon states 
because of the vanishing energy denominators; in this case the explicit diagonalization is 
mandatory. 

The mixing matrix V^a has a rather peculiar wing-like shape as illustrated in Fig. ^. 
We remind that the index a along the horizontal axis serves to numerate different quark- 
antiquark-gluon states, starting with the one with the lowest anomalous dimension at 
a = 0, and the index j3 enumerates the three-gluon eigenstates in the similar fashion. 
Dark regions in Fig. ^ correspond to large absolute values of V^a and light regions indicate 
small matrix elements. Note the complicated 'chess-board' pattern with alternating large 
and small entries. 

The most important feature that is seen in Fig. ^ is that the lowest quark eigenstate 
a = mixes significantly only with the lowest gluon eigenstate (3 = 0. In fact, we find a 
(roughly) exponential hierarchy of the matrix elements 

|Vo,o| > |Vi,o| > ... > |V[jv/2]-i,o| , (4.21) 



n r- 




a 

Figure 6: A density plot for tlie matrix |V/3q,| for N = 100. 



that is valid for all A^, see Fig. ^. Higher quark-antiquark-gluon states get mixed with 
gluons more heavily and involve many three-gluon states in an essential way. For a 
few highest quark states, a ^ N , the mixing is again simplified somewhat, but still 
involves several (in general ~ In N) gluon states, see Fig. |^. It can be shown that the 
mixing coefficient of the highest quark-antiquark gluon state ( |4.5| ) with the lowest three- 
gluon state ( [4.14| ) Vq^n vanishes in the large-iV limit and this smallness overcomes the 
enhancement due the small energy denominator in ( |4.18| ) in the same limit, cf. ( [4.17| ). 

Since the Hamiltonian ([4.1|) is hermitian, the same mixing matrix describes the mixing 
of the gluon states with the quark states. It is seen that the lowest gluon states (3 = 
0, 1, . . . get mixed with many quark states while the highest gluon states f3 ~ N/2 only 
mix with quark states with a ~ 2/3N. 

As a yet another illustration of the mixing pattern, we show in Fig. |^ the exact nu- 
merical results for the expansion coefficients of two different exact eigenstates of the 'full' 
Hamiltonian ( |4.1| ), "ifN^^Xi) and \l/iv Ar+i(xj), over the 'pure' quark and gluon eigenstates. 
The chosen states are those that reduce to the lowest quark and the lowest gluon states if 
the mixing is 'switched off'. For the lowest quark state (which is the lowest eigenstate in 
the whole spectrum) the single significant contribution comes from the lowest gluon state 
(~ 10^'^) while all other contributions are suppressed by another order of magnitude. 
For the gluon state the mixing is much larger ( < 10~^) and involves many quark states. 
This is the same structure that we observed earlier in Fig. p. Another conclusion is that 
since the mixing-induced corrections are at most of the order of 10% (for large A^ ~ 100), 
the perturbative expansion in (|4.18|) should be rather accurate and this indeed can be 
verified by the direct numerical calculation. 
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Figure 7: The mixing coefficients V/jo of tlie different three-gluon states f3 = 0,1, ... , [N/2\ — 1 
with the lowest quark-antiquark-gluon state, for four different values of A^: = 8 (triangles), 
= 20 (squares), = 50 (open circles) and A^ = 100 (full circles). 



4.3 Reduction of the mixing matrix: The two-channel DGLAP 
equation 

We have defined the transverse spin quark and gluon distributions Aqxix) and AgT{x) 
as the specific projections of the generic quark-antiquark-gluon and three-gluon opera- 
tors, Eqs. (|2.28|) and (|2.29| ), respectively. They correspond to the leading order quark 
and gluon contributions to the operator product expansion of the T-product of the two 
electromagnetic currents ( |2.1CI| ) and determine the odd moments of the structure function 
( p.l9| ) at a high scale Q"^ 



[92{n,Q')] 



2Mtw-3 



- Ag+(n, Q2) + ^Cj{n) Agrin, Q') 
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TT 



(4.22) 



Here, we took into account the relation between the gluon coefficient functions ( p.l02| ) 
and introduced a notation for the overall gluon coefficient function Cj 
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n(n + 1) 



c{n)[i; (n) +-fE+ 1]) 



(4.23) 



with c{n) given by ( p. 1041 ). It is therefore natural to exercise a model for the structure 
function in which g2{x, Q"^) can be expressed in terms of these two distributions defined at 
a low scale /i^. Such a model cannot be exact, since it is not theoretically consistent: Other 
partonic degrees of freedom are generated through the QCD evolution of Aq'j'{x, Q"^) and 
Ag^ix, Q^) down to a low scale. Assuming the absence of their contribution at two 
different scales and /i^ simultaneously is, therefore, not possible. We shall argue. 
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Figure 8: Exact numerical results for the expansion coefficients of the exact eigenstates 
of the Hamiltonian H,^ over the eigenvectors of the diagonal blocks: a) - the exact lowest 
eigenstate, (\E'7v^o|^7Vfc)' ^) ' exact lowest eigenstate in 'gluon' sector, (\E'Ar^7v+i|^7Vfc)- 
The index k enumerates the operators ordered according to their anomalous dimension; 
< A; = a < corresponds to quark operators, k = N + l + f3>N + lto gluons, 
cf. (f4.3|). The 'quark' eigenstates are shown by crosses, the 'gluon' eigenstates by open 
circles. 



however, that the admixture of the additional degrees of freedom turns out to be small 
numerically, at least for large A^. This means that the two-component model for the 
structure function g2{x, Q^) based on the parton distributions Aqt^x, Q^) and Agx^x, Q^) 
has a good numerical accuracy and can eventually be made more sophisticated (involve 
more states) when and if high accuracy experimental data become available. 

In order to unravel the qualitative structure of the evolution it is necessary to go over 
to large moments so that the number of contributing operators becomes large and 
possible systematic effects more transparent. To this end we take = 100 and, as a ffist 
step, examine the coefficients in the expansion 

and AgxiN + 3) = dx x^~^'^ AgT{^, Q'^) in contributions of the multiplicatively renor- 
malizable operators (|3.83|) . We recall that the moments Aq^{n) and Agxin) are given by 
reduced matrix elements of the local composite quark- ant iquark-gluon and three-gluon 
operators defined as 

Ag+(n) = 2{{<^l{d,) S,{z,))) , Agrin) = (($^„(9,) d,{z,))) . (4.24) 
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Replacing the nonlocal operators by the expansion over the multiplicatively renormaliz- 
able operators (p.l6|) one gets 
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Figure 9: Exact numerical results for the coefficients 



and 



N,k\^N+3/ 

b) — K^jv-i,A:l'^Ar+3)l/ll'^Jv+3ll ^^e expansion of the N = 100th moment of quark and 
gluon transverse spin distributions, respectively, in the contributions of multiplicatively 
renormalized operators. The index k enumerates the operators ordered according to their 
anomalous dimension; < A; = a < corresponds to quark operators, k = N + l + l3> 
+ 1 to gluons, cf. (|4.3| ). The 'quark' eigenstates are shown by crosses, the 'gluon' 
eigenstates by open circles. 
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a=0 



where ^ 



and \&jv-i a 'quark' and 'gluon' components of the exact eigenstate of the 
evolution kernel in the helicity basis. 

The expansion coefficients (•^'at+sI^atq) for N = 100 are shown in Fig. It is seen 
that Aq^[N + 3) is grossly dominated by the contribution of the local operator with the 
lowest anomalous dimension, which means that all other contributions to the evolution 
are small. Among the other contributions, noticeable corrections come from a few quark 
operators with the low anomalous dimensions, and a single gluon operator with the lowest 
anomalous dimension in the gluon sector. The contributions of the quark operators with 
the anomalous dimensions close the the lowest one are less interesting than the gluon 
contribution for the following two reasons: 

• Since the anomalous dimensions of important quark operators are all close to that 
of the leading quark operator, the corresponding contributions are barely distin- 
guishable by the evolution. 

• The admixture of quark operators with the next-to-the-lowest etc. anomalous di- 
mensions is roughly the same in the flavor-nonsiglet and the flavor-singlet channels. 
Because of this, one does not expect any qualitative effects. On the other hand, 
the appearance of the gluon operator is a new feature of the singlet evolution. 

Neglecting the higher quark and gluon states in ( |4.25| ) we find that for arbitrary (large) 
A^ the moments of the quark distribution receive the dominant contribution from only 



two states - the lowest quark (a = 0) and the lowest gluon (a = + 1) levels and, 
therefore, depend on only two nonperturbative parameters. As a consequence, thus 
defined moments Aq~^{n,Q'^) satisfy the two-channel evolution equation which can be 
written in the standard DGLAP form 



K(^/2/)Ag+(i/; Q') + P^^{x/y)Agr{y; Q')] . (4.26) 



In order to determine the splitting functions one can take moments so that the convolution 
in ( |4.26| ) gets replaced by a product and identify the moments of the splitting functions 
with the corresponding anomalous dimensions (with sign minus). 

Since moments of the quark distributions are obtained from the nonlocal quark- 
ant iquark-gluon operators by projecting onto the quark coefficient function ( [4.24|) , one 
can calculate the anomalous dimensions by applying $^(i9j) to the evolution equation 
( ^.14| ). Since the coefficient functions $^_,_3 and $^+3 of the quark and gluon transverse 
spin distributions Aqx^N + 3) and Aqt^N + 3) ( |2.30| ) turn out to be very close to the low- 
est eigenstates of the 'pure' quark and 'pure' gluon Hamiltonian Hq ( [4.1| ), see Eqs. ( |4.10| ), 
( [4.14| ), we may invert ( [4.24 ) and expand the nonlocal operators S^{zi) and 0^{zi) as 



iS,{z,))) = Y.^.Aq^M') + . 

n " " 



(4.27) 



where $^ and are related to the coefficient functions through the transformation 
( ^.34| ) and the dots denote the contribution of higher quark and gluon eigenstates. Then, 
substituting ( |4.27| ) into the evolution equation ( |3.14| ), one calculates the corresponding 
anomalous dimensions (for even A^) as 
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where the last factor in the second expression serves to correct for the different normal- 
ization of the quark and gluon coefficient functions. The first matrix element has been 
calculated in [p!0| , [1^ and the answer for the large- A^ expansion of E^jq is given in ( |4.12|) , 
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For the second matrix element we have 
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Figure 10: Same as in Fig. ^ but for = 10. 



Using the explicit expressions for the coefficients 0^ ^ (|3.90|) and 



Lh 
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^ (|3.97|) one gets 
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Combining together Eqs. ([4.31|) , (|3.100|) and ( p. 911) and expanding at large n 
we obtain 
-1 
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(4.31) 

A^ + 3 



(4.32) 



which agrees well with the exact expression for all n > 3. 

In order to get a closed system of equations, we have to consider the evolution of 
Agxix) as well. The coefficients of the expansion of the corresponding coefficient function 
over the basis of multiplicately renormalizable operators, Eq. ( ^.25| ), are shown in Fig. 
for = 100. We see that the gluon distribution is dominated by the contribution 
of the gluon operator with the lowest anomalous dimension in the gluon sector, and 
contributions of all other gluon operators with higher anomalous dimensions is strongly 
suppressed. Notice, however, that the contribution of the quark operator with the lowest 
anomalous dimension is small compared to the contributions of a large number ~ iV of 
other quark operators with larger anomalous dimensions. This means that although we 
can formally write 
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where Sf^^ is given in ( ^.15| ), taking into account the mixing term ~ Pj^i^x/y) is not 
justified to the accuracy that the mixing with other quark degrees of freedom is omitted. 
We suggest, therefore, to neglect the mixing of quarks to the gluon distribution altogether. 
It is this mixing that sets the limitations for the accuracy the two-component transverse 
spin model. This accuracy is in fact rather high since the sum of squares of the coefficients 
of all quark operators in Fig. ^3 is less than 2%. 

For smaller values of A^, the hierarchy of different contributions does not look so con- 



vincing, see Fig. |T0|, but qualitatively remains the same. Most importantly, contributions 
of all gluon operators other than the one with the lowest anomalous dimension remain 
negligible. 

The DGLAP equations ( [4. 261 ) and ( [4.33| ) can easily be solved going over to moments. 
Neglecting the quark mixing in ( [4.33| ) we find using ( f4.22| ) 
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""'~^P^Ij{x) and L = as{Q^)/as{fi'^). For the two lowest moments 



where 7jb(?T.) = — dx x' 
n = 3 and n = 5, expanding moments of the parton distributions at a low normalization 
scale Ag^(n, /i^) and Agxin, fi"^) over the reduced matrix elements of local operators 



^)) and {{[G]%)), Eqs. ( p.31| ) - (p.34|) , and using the explicit expressions for the 



evolutions kernels, one obtains 

^10.987/6 - 2 {{Si)) + 3 ((5°)) + 0.974 

+ L™/" [- 0.974 



(4.36) 



These results have to be compared with the exact expressions in ( ^.9] ). We see that the 
approximation advocated in this work corresponds to taking into account (for n = 5) the 
two multiplicatively renormalizable operators with a large gluon contribution (the first 
and the third lines in ( p.9|) ) and neglecting the other terms. The coefficients in front of 
the three quark local operators in the first line in ( |3.9|) are reproduced reasonably well. 
Note that values of the anomalous dimensions come out to be very close to the exact 
results. 



5 Summary and conclusions 

Based on our systematic analysis of the evolution pattern of twist-three operators we 
arrive at the following approximate two-channel evolution equation for the fiavor-singlet 
quark and gluon transverse spin distributions: 



where the flavor- singlet quark distribution is defined similar to 

Aqp^{x) = Aut{x) + Aut{x) + Adrix) + Adrix) + 



(5.1) 



(5.2) 



The Hamiltonian approach developed in this paper can be used to determine the asymp- 
totic expansion of the moments of the splitting functions at large N that translates to 
the expansion in powers of 1 — a; at large momentum fractions x. Since quark and gluon 
distributions in the nucleon appear to be decreasing strongly at large x, this contribu- 
tion to the splitting functions is the most important one numerically. Using the results 
for the moments given in the text, we derive the following expressions for the splitting 
functions: 
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Here, the first two expressions are accurate up to corrections of order 0{1 
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(5.3) 
for 
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1 and the third expression has the accuracy 0{{1 



X] 



Note that the quark 



splitting function Pqq{x) is the same in the flavor singlet and flavor- nonsinglet channels. 



cf. ( p..2|) . Remarkably enough, the obtained the twist-3 evolution kernels turn out to be 
very similar to the well-known expressions for the twist-2 DGLAP kernels. Moreover, the 
leading x — > 1 asymptotics of the diagonal kernels is the same for twist-2 and twist-3 and 
the difference occurs at the level of sub leading ~ 6{1 — x) corrections. One can argue 



following that this property is rather general and it holds to all orders of perturbation 
theory. 

To the leading logarithmic accuracy, the structure function g2{x,Q'^) is expressed 
through the quark distribution 



y 



(5.4) 



where g2^^{x, Q"^) is the Wandzura-Wilczek contribution ( |2.7|) and the gluon contribution 
arises entirely through the evolution equations (|5.1| ) (after the separation of the flavor- 
singlet part). Note that both the gluon and the quark distributions are defined by analytic 
continuation from the odd moments ra = 1, 3, . . . and must satisfy the constraints 



dx Aq^{x, 



dx Agxix, Q^) = 



(5.5) 



that follow from properties of the coefficient functions $^ in ( p.3(]| ). The Burkhardt- 
Cottingham sum rule dx g2{x,Q^) = then follows from (|5.4| ) which derivation 



involves an additional implicit assumption about the absence of subtraction constants in 
the dispersion relation for the spin-dependent Compton amplitude, see [^. 
The gluon distribution is subject to the additional constraint 



dxx^AgT{x,Q^) = 



(5.6) 



which ensures that the gluon contribution vanishes for the third moment. 

To the next-to-leading logarithmic accuracy, the twist-3 gluon contribution at large 
scales can be calculated as a finite part of the box diagram in the background gluon field 
and the result [|T^] projected onto the gluon transverse spin distribution has the form 



. + 



a. 



TT 



C^Jx/y) Agriy.Q' 



:i-x) 



X) 



1 + ln. 



X 



1 — X 



(5.7) 



- ln''(l -x) ln(l -x) -\ 

3 ^ ^ 9 ^ ^ 54 9 



+ 0{{l-x) 



where the dots stand for the leading-order quark contribution and the 0{as) quark cor- 



rections. We have argued in |jT5[ and in this paper that this projection has a high accuracy 
for all integer moments (i.e. all contributions left out by this projection are very small). 

To summarize, the set of formulas given in this Section presents a theoretically mo- 
tivated approximation for the QCD description of deep inelastic scattering from a trans- 
versely polarized nucleon. Importance of this approximation is not so much in the possi- 
bility to calculate the scale dependence, but in the identification of important transverse 
spin degrees of freedom that are preserved by QCD interaction. The formalism devel- 
oped in this paper is very general and can be applied to the study of other higher-twist 
distributions. It is based on the operator product expansion and conformal symmetry of 
the QCD Lagrangian. This technique turns out to be very effective in dealing with the 
operator renormalization for large and in many cases a WKB-type expansion in 
can be constructed analytically. Going over from the moments to the momentum fraction 
space involves analytic continuation and in general may be quite complicated, see for 
a discussion. Because of this, the small— x behavior of twist-3 (and higher-twist) parton 
distributions presents a nontrivial problem and deserves further study. 
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Appendices 



A A 



A Appendix: Calculation of the anomalous dimen- 
sions 



In this Appendix we describe the calculation of the energy of the lowest quark and gluon 
levels, Eqs. ( [4.11|) and ( [4.16|) , respectively, and the energy of the highest quark level. 



Eq. (O). 

The energy of the lowest gluon level is defined by Eq. ( 4.15|) and involves the ex- 



pectation value of the Hamiltonian H^h over the gluon state and its norm. To 

calculate the both, we shall use the expansion of the wave function over the confor- 
mal basis. An arbitrary three-gluon state can be characterized by total conformal spin 
J = n + 1/2 = N + 7/2 and the conformal spin in a certain two-particle channel. De- 
pending on the particular choice of the two-particle channel one can define three different 
sets of basis functions Y^^l^^^^j^^{xi), ^i+i/2,fc+3(a;j) and ^„+\/\fc+3(a;i)- The functions be- 
longing to each basis are linear independent whereas the ones belonging to two different 
sets are related to each other through a linear transformation 

^ n+in.k+Zy-^-^'l ~ / .^^km ^ n+ll2.m+Z\-^^) ~ / \ ^) "fcm n+l/2,m+3 '^'^ V ' \^-^) 
m m 

with = (>;+\/\fc+3l>;+?/2,™+3) being the Racah 6j-symbols. 

The gluon coefficient function in the helicity representation can be expanded in any 
one of the three conformal basis so that one gets three sets of the expansion coefficients: 



n— 4 n— 4 n— 4 

E, (31)2 T^(31)2 / \ _ ,(12)3 ^{12)3 / n _ , (23)1 ^^(23)1 

'rn,k ^n+l/2,k+3\^i) ~ 2-^'rn,k n+l/2,k+2,\^i-) ~ 2-^Vn,k n+l/2,A;+3' 
fc=0 k=0 k=Q 



These coefficients can be calculated similar to ( |3.96| ) and are given by 

{n + 2 + k){n + A + k'^/2 + 5 k/2) 



,(12)3 i,(31)2 , 

V =-<fc =Kk 



{k + l){n-2-k) 



(23)1 
n,k 



Xi). 

(A.2) 
(A.3) 



for even k, and 

, (12)3 _ ,(31)2 _ 1^(23)1 {k + A){n + 3 + k) 

Kk -Kk --l^Kk - hn,k - 3 - k) ^^-^^ 

for odd k. Here, the normalization factor hn^k is defined as 

2 _ {k + 1) (-n + 3 + k) {-n + 2 + k){2k + 5) T\n) 



The relations between the coefficients respect the symmetry properties of The 

f 23) 1 

coefficients (j))^ ^ vanish for even k to ensure the antisymmetry of the state under the 
interchange of gluons with the same helicity, X2 ^ x^. The sum of the coefficients in the 
three different channels vanishes for arbitrary k since the three-gluon state is annihilated 
by the operator 1 + P + P"^ = [P^ — 1)/(P — 1) with P being the operator of cyclic 
permutations. 

The 'spherical harmonics' ^^^^^^(a^i) form an orthonormal basis on the space of the 
coefficient functions endowed with the scalar product ( p.39|) , {Y^^^^'^\Yn'^m'^) = 6k^m- The 



same is true for the states F 



(12)3, 



and Y" 



(23)1 



. Using this property we may calculate 



the norm of the gluon state in one of the three equivalent forms 

n—4 ^ n—4 „ n—4 



1$ 



h\\2 



/,(12)3 
^n,fc 



/,(23)1 
^n,fc 



/,(31)2 
^n,fc 



1 2 



(A.6) 

fc=0 k=0 k=0 

Substituting the expressions for the expansion coefficients ( [A.3| ) and ( |A.4| ), one gets a 
finite sum over two-particle conformal spin k whose evaluation leads to Eq. (|3.100|) . 
Expanding the result at large n one obtains 



1$ 



hll2 



r^(n) 



1 + -(2 Inn + 

n 



(A.7) 



1 [2nj 16 

where we factored out the ratio of the F— functions for later convenience. 

In order to calculate the expectation value of the evolution kernel defined in ( p.74| ), we 
notice that Hhh has a two-particle structure and can be split in three contributions each 
of which only depends on the operator of the conformal spin Jab in a given two-particle 
channel 

{^'jH,,\^t) = Nc {^tm,,iJ,2) - K„(Jl2)|0 

+ Ar,«|2[/,,(J3i)-K„(J3i)|0 
+ 2N,{^':,\U,,iJ,s)\K) -b\\K\\' ■ 



Since the operator J12 is diagonal in the conformal basis Y 



(12)3 
n,fc 



Xi 



(A.8) 

(by construction), it 

is natural to evaluate the first term in ( A.8| ) by expanding the gluon state over this 
particular basis. By the same token, the second and the third terms in ( [A.8| ) are most 
easily calculated using the expansion over Y^^^^'^' 
way one arrives at 

n—4 

= 2iv,j][/,,(A:+3) 



Xi 



and Y^^^^^{xi)^ respectively. In this 



"IT? 



+ 



+ 



k=0 

n-A 



-^c5^K„(A; + 3) 



611$ 



h\\2 



(A.9) 



k=0 



where we have used that the conformal spin in the two-gluon channel is equal to 2jg + k = 
A; + 3. Using explicit expressions for the kernels Ugg and Vgg defined in ( ^.77] ) and the 
expansion coefficients ( |A.3| ) and ( |A.4| ), one can rewrite ( |A.9| ) as a finite sum over the 
two-particle spin k. The sum involving Vgg can be calculated analytically while the sum 
involving Ugg can be expanded in inverse powers of 1/n at large n. After some algebra 
one arrives at 



n—4 



■'n,k 



+ [0n,A 



k=0 



r^(n) 
r(2n)T 



3 
2 



n 



Inn 



and 

n-4 



fc=0 



E Ugg{k + ?>)(^ 



Inn 7s + 



TT 

12 



2 ^ n 



2(lnn + 7s)' - -(Inn + 1e)]+0 



(A.IO) 



(A.ll) 



In n 



n^ 



Substituting these expressions into ( [A.9| ) and combining them with ( |A.7| ) one arrives at 
the expression for the energy of the lowest gluon state given in the text, Eq ( [4.15|) . 

Calculation of the energy of the lowest quark level, Eq. ( [4.11| ), goes along the same 
lines. The quark coefficient function $^ is defined in Eq. ( p.22|) on the hyperplane X1+X2+ 
X3 = corresponding to the kinematics of the forward scattering, and can be continued 
to arbitrary values of the momentum fractions Xi using the conformal symmetry. Its 
expansion over the conformal 'spherical harmonics' in the three different two-particle 
channels looks like 



n— 3 



n— 3 



n— 3 



$9 (^.) _ V (.(12)3^(12)3 (x.)-y" ^(23)ly(23)l 



,k+5/2 



\P,„(31)2^(31)2 



k=0 



k=0 



k=0 



where e.g. Y 

c5 n 



(12)3 

+l/2,fc+5/2 



Xi 



is defined by the general expression ( p. 611 ) with ji 



(A.12) 
33 = 1, 



J2 = 3/2, the total conformal spin equal toJ = ?T, = A^ + 3 and the two-particle conformal 
spin in the (12) channel given by ji2 = ji + j2 + k = 3 + k. 

The explicit expressions for the expansion coefficients can be obtained using the rep- 
resentation similar to (|3.89| ). The result is: 



(12)3 ^ 2(-l)'' ^ + ^ + ^) 



(23)1 
Vn,k 

(31)2 
Vn,k 



<in,k 



k + l 

-2 (-l)"(fc + 3) (n + 2 + fc)g„,fc 
- (3 + 2 A;) (fc + 2) 



(A.13) 



X 



\n—k 



ln + 2 + k (-1)" + ln + l + k 



n-2-k 



+ 



n — 1 — k 



with the normalization factors 



^n,k 



{k + l){n-2-k) 



-2, 

Qn,k 



4:{n + 2 + k){k + 2){k + 3) T{2n) 
^ (n-1 



k){k + 3) 2 

:Qn,k 



(A. 14) 



{3 + 2k){n + l + k) 

Using the expansion in ( |A.12| ) one finds three equivalent representations for the norm of 
the quark state 

n— 3 n— 3 n— 3 

k=0 k=0 k=0 

Substituting explicit expressions for the expansion coefficients and performing the sum- 
mation one arrives at Eq. ( p.91| ). The expansion of the norm at large n reads 



|$9||2 



n 



T{2n) 4 



1 + —(1 - 4 Inn - + 0{l/n^ 



(A.16) 



The expectation value of the diagonal quark evolution kernel Hqq defined in ( 3.69|) can 
be written in the form analogous to ( |A.8| ) 



A i-r 



+ ($^|iv,f/W(j23) - —u!^i\j,,)\^i) 

Expanding the quark state ( |A.12| ) over the suitable conformal basis one obtains 

n-3 / 2 

{K\H,,\K) = Et^iT]' + 5/2) - -rrViPik + 5/2) 



(A.17) 



k=0 
n— 3 



+ + 5/2) - -^f/«(A; + 5/2)) (A.18) 



where we have replaced the operators of the two-particle conformal spins by the corre- 
sponding eigenvalues Jab = ja + jb + (n — 3). The necessary expansion coefficients and 
explicit expressions for the evolution kernels entering this expression are given in ( A.13| ) 
and ( |3.73| ), respectively. 

It turns out that the part of the sum in (|A.18 ) proportional to Nc can be calculated 
exactly 



n— 3 



E[^iT]'n?n^+5/2) + bS)?t/S)(fc+5/2) = mrUn) + i - i + 27^ ). (A.19) 

k=0 ^ 



The 1/Nc correction to the sum (|A. 18|) is given by a ratio of rather complicated sums 
and can easily be expanded at large n leading to 



n— 3 



g ||2 



k=0 
n-3 



4 6 



C(lnn/n^ 



(A.20) 



A:=0 
n—3 



Y.l^'nTMi'i^ + 2) = ||<|.^f (Inn + 7,. - 1 + 0{\n' n/n')) . 



k=0 



Finally, the n/— dependent contribution to the ( |A.18| ) is given for odd n by 



(31)2n2 



['Pnfi 



(n- l)(n + 2) T^{n) 
(n + l)(n-2)r(2n) 



m'{0 + O{l/n')) 



(A.21) 



Combining ( [A.19| ), ( [A.20| ), ( [A.21| ) and ( |A.17| ) we obtain the expression for the energy of 
the lowest quark level given in Eq. ( |4.11| ). 

The eigenstate of the highest quark level can be approximated by the expression in 
(|]6|). Calculating the corresponding energy (^^ ) and using ( p.69| ), one gets 



(A.22) 



A n 



Using the explicit expression ( p.7(]| ) for the Hamiltonian Hs+ and taking into account that 
^iVfc=o symmetric with respect to the interchange of the quarks, (P13 — 1)Y^^^q = 0, 
one obtains 



^In = 0% - —Uil\Ji3 = 2) + (YloiyiJuWl,) , (A.23) 



2 

where the notation was introduced for the hnear combination of the kernels (|3.73|) 



ViJ) = {V^fiJ) + Uj^fiJ)) - |- KpiJ) + U^l\J)). {AM) 

In order to calculate the matrix element entering ( |A.23| ) we use the Racah decomposition 
to expand over the conformal basis in the quark-gluon channel 

YU^^) = J2 ^2,n^+5/2(A^ + 7/2) Y^'ll),^^^^^,{x,) . (A.25) 

m 

Taking into account that Uql\2) = one finds 

2 ^ 

^NN = i;nf + Y, + 5/2) [^^2,m+5/2(A^ + 7/2)]^ . (A.26) 

m=0 

Finally, using the explicit expressions for the Racah symbols ( p.3|) and the two-particle 
evolution kernels ( |3.73| ) and and performing the summation in ( |A.26|) one arrives at the 
result given in Eq. 



B Appendix: Racah symbols for the SL{2, R) group 

In this Appendix we derive the explicit expression for the Racah 6j— symbols Qjj' defined 
in (Ipl) 

yif"i^^)= E (B.i) 

ji+j2<j'<J-j:i 

where the basis functions Yjj are given in Eq. ( |3.61| ). In order to find the coefficients 
Qjj/{J) it is sufficient to set the three variables Xi in the Eq. ( |B.1| ) to the following values: 
X\ = — X, X2 — X and X3 = 1. Using the explicit expressions for the basis functions ( |3.61| ) 
one gets for 



( i)J-.i-h^, r(j+j3-ji)r(j-j+j2) . . 

^ "^^r(2j3)r(2j,)(j-ji-j3)!(J-j-j2)! ^ ^ 

X 2F1 (ji + J3 - J, 1 +33 -31 -3, 2j3 ; X) 2F1 {-J +J +J2 , J +j +32 - 1 , 2j2 ; x) 

^ V o ( n "^^^ (-^)^"^'^"^'^ nJ+3'-33) 

.'4^. (^'-^■i-^2)!(J-j'-j3)!(2j'-l)r(j'+ii+j2-l)- 

The product of the two hypergeometric functions in the l.h.s. of this indentity is nothing 
else but the generating function for the Racah polynomials (see for details). Expand- 
ing the hypergeometric functions and comparing the terms in the l.h.s. and the r.h.s. 



A r\ 



of ( |B.2| ) with the same power of x one derives after some algebra the following explicit 
expression: 



1 



r(J-ji-j2-j3+i) 



f{J, f, ii, j2, is) f{J, j, 31,33,32) 



r(2ji)r(j+ji+j2+j3-i 

-f +31+32, j'+jl+j2-l, -j+jl+j3, j+jl+j3-l 

2ji, -J+31+32+33, J+3i+32+33-'^ 



(B.3) 



where 

f{J, 3 ,31,32,33) 



(2/-1 



r(/+ji-j2)r(/+ii+j2-i)r(j-/+j3)r(J+/+j3-i; 



r(j'+j2-ji)r(j'-ji-j2+i)r(j-j'-j3+i)r(J+j'-j3) 



1/2 



(B.4) 



C Appendix: Evolution kernels in the conformal ba- 
sis 

To solve the Schrodinger equation (|3.58|) we expand the eigenstates over the conformal 
basis in the quark-antiquark-gluon and three-gluon sectors, Eq. (|3.66|) . In this represen- 



tation the evolution kernels are given by real and symmetric matrices acting on the vector 
of the expansion coefficients (m^^,, u^n k) defined in (p.66|) . In this appendix we work out 



the explicit form of these matrices and discuss their properties. 

The diagonal quark evolution kernel ( p.69| ) is described in the conformal basis Y^^^ 
by a square matrix of dimension ^q = N + 1 

[H,,U = {Yl,\HM,J = {Yl,\Hs.\YlJ + -^6,,^6,,o (C.l) 

with < k,m < N. In turn, the Hamiltonian Hs+ is given by the sum of pair-wise 
kernels ( p.73[ ) depending on the conformal spins in the different two-particle channels. 
We recall that the operators J31 are diagonal by definition on the space of the states ^ 
so that 

{Yl,\V{Jsi)\YlJ = V{k + 2)6,^ (C.2) 

for an arbitrary V. On the other hand, contributions to the Hamiltonian that are func- 
tions of J12 can be expanded over the conformal basis in the two-particle (12)— channel 

N 

{YN,k\^{Jl2)\YN,,J =^y{l + ^/'^){YN,k\YNll)2,l+5/2)(^N+y2,l+5/2\^N,^^ • (^-3) 

The similar representation exists for l^(J23)- The scalar product of the ^-functions 
belonging to different conformal basis is given by the Racah 6 j— symbols defined in ( |B.1| ) 
and ( |B.3| ). In particular 



(4+7/2,+5/2l>^^,™) = ^m+2,l+,/2{N + 7/2) = Q^^] , (C.4) 
(^1+7)2,^+5/2 1 = (-l)'+"fir„+2,/+5/2(iV + 7/2) = (-1)'+™^^! . (C.5) 



Taking into account these relations one finds 



N 



(9)0(9) 

kl ; 



1=0 
N 



(g)n(9) 

kl 



(C.6) 



(C.7) 



1=0 



Using these identities and the exphcit expression for Hs+ given in Sect. 3.5.1 one finds 
the following matrix representation of the diagonal quark kernel 



qq\ km 



N 



4,m + ^ ^ 



(g) 0(5) 

kl 



1=0 



X ( iV,VW(/ + 5/2) - —Vip{l + 5/2) + (-1)^+™ 



iV,f/(;)(/ + 5/2)--f/(J)(/ + 5/2) 



The explicit expressions for the Racah symbols and the evolution kernels are given in 
Eqs. ( |B.3D and ( ^.73| ), respectively. 

The gluon diagonal kernel in the helicity representation ( |3.74D is given in the conformal 
basis Y^_^ ^ by a square matrix of dimension ig = N/2 



hh\ km 



(C.9) 



with < /c, m < [iV/2] — 1. To obtain the explicit expression one has to follow the 
steps similar to that in the quark case. Introducing the notation for the relevant Racah 
symbols 

{yS'^7)2,JyN-l,m) = fi™+3,/+3(iV + 7/2) ^ Q^^} (C.IO) 

and using the explicit expressions for the gluon kernels ( |3.75| ) one finds 



[Hhhlkm = Sk,m[N, {2Ugg{k + 3) - Vgg{k + 3 ) ) " 6] 



(C.ll) 



[N/2] 



+ NcY. 4?^[2^..(^ + 3)(1 + (-1)'= + "^) - Vgg{l + 3)] 



1=0 



The explicit expressions for the Racah symbols and the functions Ugg, Vgg are given in 
Eqs. (|B.3|) and ( p. 771) , respectively. 

Finally, the explicit expressions for the off-diagonal kernels Hgh and H^g are given in 
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